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ABSTRACT 

Transitional disks are protoplanetary disks around young stars, with inner 
holes or gaps which are surrounded by optically thick outer, and often inner, 
disks. Here we present observations of 62 new transitional disks in the Orion 
A star-forming region. These were identified using the Spitzer Space Telescope's 
Infrared Spectrograph and followed up with determinations of stellar and ac- 
cretion parameters using the Infrared Telescope Facility's SpeX. We combine 
these new observations with our previous results on transitional disks in Taurus, 
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Chamaeleon I, Ophiuchus and Perseus, and with archival X-ray observations. 
This produces a sample of 105 transitional disks of "cluster" age 3 Myr or less, 
by far the largest hitherto assembled. We use this sample to search for trends 
between the radial structure in the disks and many other system properties, in 
order to place constraints on the possible origins of transitional disks. We see a 
clear progression of host star accretion rate and the different disk morphologies. 
We confirm that transitional disks with complete central clearings have median 
accretion rates an order of magnitude smaller than radially continuous disks of 
the same population. Pre-transitional disks — those objects with gaps that sep- 
arate inner and outer disks — have median accretion rates intermediate between 
the two. Our results from the search for statistically significant trends, especially 
related to M, strongly support that in both cases the gaps are far more likely to 
be due to the gravitational influence of Jovian planets or brown dwarfs orbiting 
within the gaps, than to any of the photoevaporative, turbulent or grain-growth 
processes that can lead to disk dissipation. We also find that the fraction of Class 
II YSOs which are transitional disks is large, 0.1-0.2, especially in the youngest 
associations. 

Subject headings: accretion, accretion disk — planetary systems: protoplanetary 
disks — stars: pre-main sequence — infrared: stars — X-rays: stars. 



1. INTRODUCTION 

Transitional disks (TDs) are protoplanetary disks around young stars which are optically 
thick and gas-rich, but which have AU-scale radial gaps or central clearings in their dust 
distribution. Usually the gap, which is depleted of small dust grains, is revealed as a deficit 
in the host T Tau star's infrared excess, with respect to its siblings. Such disks are thought to 
represent an evolut ionary stage between Class II and Class III young stellar objects (YSOs)E 



Strom et al.l (119891 ) first drew attention to these objects, using examples with small — almost 
photospheric — infrared excess at short mid- infrared (5 — 50//m) wavelengths but switching 
sharply to large and strong excess at longer wavelengths. 



As in lEvans et al. ( 20091 ) . we mean by Class II YSO an object with 2-20 /im spectral index between -1.6 



and -0.3, which physically corresponds to a (T Tauri) star surrounded by a (radially-continuous) accretion 
disk which is not viewed close to edge-on. Similarly, by Class III YSO we mean an object with 2-20 
/im spectral index less than -1.6, corresponding to a pre-main-sequence star with little or no circumstellar 
material. 
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During the past decade, there have been significant improvements on the s tudy of YSQs 



and protoplanetary disks based on data from the Spitzer Space Telescope (I Werner et al. 

i 



2004j ) launched in 2003. Detailed studies of disk structur e have been possib le with mid-IR 
spectra taken using from the Infrared Spectrograph (IRS; iHouck et al.l 12004 ). In particular, 



there has been major progress revealing the detailed disk structures and uncovering a variety 
of transitional disks. 



After the distinctive spectrum of CoKu Tau/4 was discovered by iForrest et al.l ( 12004 ). 
several transitional disks in th e Taurus star-forming region were stu died in deta i l with 
self-consistent disk modeling b vlD'Alessio et al.l (120051) (CoKu Tau/4). ICalvet et al.l (120051 ) 
(DM Tau and GM Aur ) and lEspaillat et al.l (12007 ajJbl. 1io08h (LkCa 15 and UX Tau A). 
Calvet et al.l tood ) and iD'Alessio et al.l fepQsh suggested that DM Tau and CoKu Tau/4, 
respectively, can be explained with an empty central cavity in the disk surrounding the cen- 
tral star, but GM Aur requires a disk structure with a gap separatin g an inner opticall y 
thin disk and an outer opti cally thick disk to match the observed SED (ICalvet et al.l 120051 ). 
Espaillat et al.l (120081 120101 ) confirmed that LkCa 15 and UX Tau A have gaps between opti- 
cally thick inner and outer disks. We will refer in the following to these three types of disks 
as classical transitional disks (CTDs; e.g. DM Tau), weak-excess transitional disks (WTDs; 
e.g., GM Aur) and pre-transitional disks (PTDs; e.g., LkCa 15), respectively. 

There is generally good agreement with the gap sizes inferred from SED modelling, and 
direct observations of the outer edges of the gap s. Gaps in some 15 transitional disk s have 
been confirmed by su 



Andrews et al. 2009 



jmilli meter interferometry (IPietu et al.ll2006l ; iHughes et al.ll2007l . 12009 



20111). and one of t hese, LkCa 15, also in near-infrared light scattered 



by the outer disk ( iThalmann et al.ll2010l ). The resolved inner di sk of a pre-transitional disk, 
T Cha, also confirms its gapped structure (jOlofsson et al.l 120 111 ). 



The definition and selection criteria of transitional disks have been far from homogeneous 
considering all objects called 'transition' or 'transitional' objects from other works. Many 
authors define transitional disks rather loosely, with diverse nomenclatures^]. Circumstellar 
disks with a deficit of dust emission from the inner disk and large 30/1 3 nm flux ratio indi- 
cati ng an optical l y thic k outer disk have been designated "cold disks" by lBrown et al.l (120071 ) 
and lMerm et al.l (120101 ) . A more relaxed definition includes disks with a much smaller excess 
for wavelengths beyond > 13 fim representing a depleted or settled outer disk compared t o 
our definition of transitional disks. Those disks are named "anemic disks" JLada et alil2006h 
"homologously depleted disks" (jCurrie et al.ll2009l ) , or "weak-excess disks'o (jMuzerolle et al. 



2 The definition of diverse nomenclatures are summarized in Evans et al.l (|2009) 

3 The definition of WTD (weak excess transitional disk) in this work and in iMuzerolle et"al1 |201ol ) are 
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2010). Such S 



Espaillat et al. 



Ds c an result from dust settling to the disk midplane ( iLuhman et al.l 12010 



2012). or outwar d truncation of a disk by gravitational interaction with a 



companion (IMcClure et al.ll2008l ). Therefore, it is important to clarify by which definition 
and selection criteria the sample of transitional/transition disks are selected, especially if the 
sample is to be used for searching for any trends to understand the properties of transitional 
disks and their origin. 

Several physi cal mechanisms have been suggested to ex p lain inside-out disk dispersal: 



photo evaporation (jClarke et al.ll2001l ; lAlexander et al.l 2006a. r ilGorti et al. 



2012 ); grain growth/coagulation/settling ( Dullemond Sz Dominiki 



2005 



20091: 



Owen et al. 



CieslJ2007f): inside- 



out c l earing by Magneto- Rotationa l Instability (MRI) ( IChiang &: Murray-Clayil2007tlSuzuki et al. 

201o |; |Perez-Becker fc Chia^3l201lh: dynamical effects by stellar /substellar companions (jLubow &: D'Angeh 
2006]); giant planet f o rmation dMarsh fe Mahonevlll992l: iRice et all [20031 buillen et al.ll2004j ; 



Varniere et ahlbood : iLubow & D'Angelol bood : Izhu et al.ll2~01~l[ boi2h . To distinguish and 



understand which mechanisms are dominant for the origin of transitional disks, not only 
flux density information at all wavelengths for a specific target but also empirical trends 
from observational data in a large sample are indispensable. 

There have been several efforts to use statistical trends among several properties of 
transitiona l disks and their host stars to constrain models of the origins of these objects. 
Kim et al.l (120091 ) analyzed IRS spectra of 13 transitional disks in Taurus and Chamaeleon 
(1-3 Myr) and fou nd a strong correlation between stellar mass and outer gap radius, R wa u- 
Merin et al.l ( 120101 ) collected broadband SEDs of some fifteen "cold" disks from several young 



association s. Including some objects reported in the literature (e.g. iBrown et al.l 12007 



Kim et al.l |2009|) to improve the sample size, they found the gap radius to scale linearly 
with M*, and to be significantly correlated with disk mass, with more massive disks tend- 
ing to have larger holes. In the somewhat older (4-12 Myr) clusters, Tr 37 and NGC 7160, 
Sicilia-Aguilar et al.l ( 120101 ) found the median accretion rates in the transitional-disk systems 
to be about an order of magnitude smaller than in Tau rus and Chamaeleo n, but similar to 
these associations' normal Class II YSOs, in contrast to lNajita et al.l (120071 ) who found that 
accretion rates for transitional disks in the Taurus-Aurigae association are systematically a 
factor of 10 smaller than the normal Class lis. 

In this paper, we present IRS spectra of 62 newly selected transitional disks in the 
Orion A star-forming region and near-IR spectra of 52 of them from SpeX/IRTF followup 
observation. Orion A consists of the Lynds 1640 dark cloud which includes the Orion Nebula 
Cluster (henceforth ONC) and the Lynds 1641 (henceforth L1641). This cloud stretches over 



different. 
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~ 30 deg 2 on the sky with Declination from — 4 d 30 m to — 9 d . In this paper, we assume objects 
to the north of <5j2ooo = ~^> d belong to ONC while objects to the south belong to L1641. 
The distan ce to the Orio n A co mplex has been estimated to be somewhere between 300 
to 600 pc (IMuench et al.l (120081 ) and refe rences therein ) . He re, we adopt the distance to 
Orion A as 414 pc based on the st udy by iMenten et al.l (120071 ). We take the median age of 
ONC members as less tha n 1 Myr ( Hillenbrand 1 199 7a ) and that of LI 641 members as 1 Myr 
jGalfalk fc OlofssoiJl2008h . 



We consider 105 transitional dis ks not only of Orio n A including ONC and L1641, 
b ut also of Taurus (henceforth Tau; iFurlan et al.l 120111). Chameleon I (henceforth Cha 



I; iManoi et all 1201 



), Ophiuchus (hen ceforth Oph: iMcClure et al.l lioioh . and NGC 1333 



Arnold et al.ll2012l ). From this very large sample of transitional disks 



(henceforth N1333; 

selected homogeneously from disk properties measured with IRS spectra, and comprehen- 
sive information on mass accretion rates from SpeX spectra and X-ray luminosity, we are 
now able to search for trends between stellar parameters and disk parameters with robust 
statistics. We anticipate this work to help in understanding mechanisms responsible for the 
origins of transitional disks. 

In Section 2, we present IRS and SpeX observations and data reduction procedures. 
In Section 3, we describe the measurements of stellar properties and mass accretion rates. 
We describe the extraction of the outer gap radius, R wa u, after explanation of how we 
selected transitional disks and their sub-classification in Section 4. We discuss the fraction 
of transitional disks and age trends in Section 5. In Section 6, we explore correlations and 
trends of the transitional disks properties and compare them to those of radially-continuous 
disks where possible. Then, we review and discuss what these findings mean and how 
they can be used to help our understanding of the origin of transitional disks in Section 7. 
Summary and conclusions drawn from our finding follows in Section 8. 



2. OBSERVATIONS AND DATA REDUCTION 

We observed the Orion A star-forming region between 2006 November and 2007 October, 
during Spitzer-IRS campaigns 36, 39, 40 and 44. In all, our targets were 555 objects selected 
on the basis of their Spitzer-IRAC and MIPS colors to belong to YSO Classes 0, I, flat- 
spectrum and II, wi th flux densities in excess of 2 mJy at 8 and 15 mJy at 24 /im 



( IMegeath et al.ll2012l ). Of these, 303 objects were classified as Class II based on their spectral 
index between 4.5 /im (IRAC channel 2) and 24 /im (MIPS channel 1) or have the colors 
of transitional disks seen in other regions; they do not show evidence of envelopes in their 
IRS spectra, and thus their infrared excess is due to circumstellar disks.. We observed 241 
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objects (114 in L1641; 127 in ONC) with the IRS with full IRS wavelength coverage of 
5-37 /im in the low resolution mode. We observed 62 additional objects located close to 
the Trapezium region with only partial wavelength coverage (5-14 /im), as the detectors 
for the longer wavelengths would have been saturated by the bright background emission. 
We estimate that our sample is complete for star-and-disk dominated objects for which the 
host-star spectr al type is M4 or earlier. Analysis of the full sample will be presented in 
more detail by iKim et al.l (120131 in preparation). The complementary samp le of objects 
with e nvelopes - Class 0, I and flat-spectrum objects - will be discussed by iPoteet et al. 
(120131 . in preparetion) . 



2.1. IRS /Spitzer 

The IRS spectra of the 62 transitional disks were selected from among the 241 objects 
which were observed with the IRS low resolution modules (A/AA ~ 90; Short-Low (SL):5.3- 
14 /im; Long-Low (LL): 14-38 /im), based on the selection criteria described in Section 4. 
Each object was observed at two nod positions separated by one-third the length of the slit 
in the staring mode. To extract the spectra, we used version S15.3 of the basic calibrated 
data (BCD) product from the Spitzer Science Center IRS pipeline for both SL and LL. 

To fix bad, hot and rogue detector-array pixels before extracting objects from the 2D 
spectral images, we generated a set of "grand rogue masks" for Orion A data. A grand 
rogue mask for the general data reduction process is generated by combining rogue mask 
files available up to a recent IRS campaign supplied from the Spitzer Science Center, and 
adding additionally identified rogue pixels from data images. However, we applied additional 
special treatment to make LL grand rogue masks for these data because (1) the LL array 
exposure to cosmic rays had increased continuously up to campaign 44 ^, which caused 
the S/N to decrease and the number of rogue pixels and hot pixels to increase, (2) fluxes 
from our objects are fainter than similar sources in other nearby star-forming regions due 
to the greater distance to OriA, (3) the large number of rogue pixels and hot pixels severely 
degraded the signal-to-noise ratio of our targets, and introduced spectral artifacts when 
many rogue pixels were grouped in clusters. To avoid rejection of too many pixels which 
might be perfectly good, we chose only the rogue pixels which appear in the campaign rogue 
masks more than 10 times through campaign 44. To these we added hot pixels selected by 
eye from LL 2D images from campaign 39 and campaign 44. The usage of grand mask files 



4 In IRS campaign 45 the bias voltage on the Long Low array was reduced from 1.8 to 1.6 volts and the 
array temperature was reduced from 4.4 to 4.1 K) 
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including rogue pixels and hot pixels cannot perfectly clean bad pixels. However, we tested 
several version of grand rogue masks for our LL data and chose a set of grand rogue masks to 
apply for all of OriA ClassII data. We fixed those rogue pixels in the grand rogue mask (and 
the permanently bad pixels) by interpolation in the spectral direction from nearest-neighbor 
good pixels. 

Spectral extraction of Class II YSOs in Orion A is generally much more challenging 
than in other nearby star-forming regions. This is due not only to the source faintness but 
also to the range of spatial structure in sky emission from the Orion HII regions, and the 
high stellar density. We used four source extraction methods to derive the final point-source 
spectra, choosing among them to optimize the rejection of emission from the sky and other 
nearby point sources. 

As a first, basic source extraction step, we used an automated extractor ("auto") based 
upon the IRS instr ument team's Spectral Modeling, Analysis and Reduction Tool (SMART; 



Higdon et al.ll2004j ). In auto we extracted sources from the uniformly- weighted signal along a 
tapered column, 3-5 pixels wide. In this step we used two versions of background-subtracted 
images: one prepared by subtraction of the two nods ("off-nod"), and the other by subtrac- 
tion of the sky spectrum in each grating order obtained while the target was being observed 
in the other order ("off-order"). If there were no serious sky background issues or no ad- 
ditional sources in the 2D images, the spectra taken from images subtracted by off-order 
or off-nods backgrounds are very similar. We examined each sky-subtracted image and the 
resulting spectra. If large sky-subtraction artifacts (e.g. spectral lines in the sky, contami- 
nation from nearby point sources) remained in both auto products, we re-reduced the data 
in SMART by using the same tapered column extraction as before, but with sky removed 
by fitting and subtracting a Oth- or 1st order polynomial to emission outside the span of the 
target. We designate the tapered column extraction with subtraction of a polynomial sky 
selected manually as "man". 

Usually One of three tapered column extractions gave artifact-free results even when 
the sky emission had complicated structure, but issues remained for targets with a neighbor 
closer than four pixels away along the slit. In those cas es, we resorted to optima l point-source 



extraction, using both the AdOpt script in SMART (ILebouteiller et al.ll2010l ) and another 



routine, OPSE, developed by our group ( Tayrien fc Forrest] 20131 . in preparation). They are 



complementary: AdOpt employs an empirical point response function (PRF) and can fit 
multiple objects along the slit; OPSE uses an analytical PRF and can make corrections to 
the spectrum for pointing errors, and thus was useful for extracting our target sources, as it 
can account other sources in a image. 



For the flux calibration of the spectra, we multiplied relative spectral response functions 
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(RSRFs), which were generated by dividing a template spectrum (M. Cohen, private com- 
munication; J. Marshall, private communication) by a calibrator's spectrum, extracted in 
the same way as the target spectrum. Our photometric standards were a Lac (A1V) for SL, 
£ Dra (K2III) for LL, and Mrk 231 (assumed to radiate as a power law) for LL wavelength 
greater than 32 /iin. 

We compared the spectra from all different versions of source extraction for each object, 
then we selected the final spectrum based upon freedom from artifacts. During the process, 
we combined spectra from different versions of methods if necessary to get the best spectra. 
For example, we use the combined spectrum for OriA-19: SL from AdOpt and LL from 
OPSE. Reduction choices for the final selection of spectra are noted in Table [TJ 



2.2. SpeX/IRTF 



52 out of 62 transitional objects in Orion A were observ ed at Near-IR (0.8 -2.4 /xm) 
wavelengths with the medium resolution spectrograph SpeX ( IRayner et al.ll2003l ). on the 
NASA Infrared Telescope Facility (IRTF) on Mauna Kea during the 2010A, 201 1A , and 
201 IB semesters. 

We observed our Orion A transitional disks with the Short-wavelength Cross-dispersed 
mode (SXD) on the nights of 15-17 Feb. 2010 and 25-27 Feb. and 6-10 Nov. 2011. We 
obtained spectra with various slit widths of 0.3", 0.5" and 0.8" for observations during the 
2010A semester depending on the seeing conditions of each night. We used only the 0.8" 
slit width for the observations 2011 February because the weather and seeing were generally 
poor. Among the objects observed during the 2010A semester, 8 lacked reliable spectral 
types in the literature, so we used their spectra for the determination of spectral type as 
well as accretion rate. Of the 8 objects, OriA-34 was only observed with the 0.3" slit to 
get the highest resolution (R = 2000) available at SpeX in order to determine the spectral 
type as well as measure the mass accretion rate under good conditions (seeing ~ 0.3"). The 
spectra of the other 7 objects were obtained with the 0.8" slit (R ~ 800-1200), sufficient 
to narrow the spectral type range down to 2-3 sub-types by comparin g absorption features 
of the Na I, Al I, M g I, Ca I and CO to those of template spectra (IRayner et al.l l2009al ; 
Cushing et al.ll2005l ). The details on spectral type determination are described in §2.2.1. For 
targets with very close neighbors, we oriented the slit so as to observe both simultaneously. 
For airmass greater than 2 we kept the slit orientation fixed to the parallactic angle. The 
SpeX/IR TF observation log is given in Table |2j We reduced our spectra with the Spextool 
package (Cushing et al.l 120041 ). and the flux calibration and telluric absorption correction 



( IVacca et al.l 120031 ) were done with a spectrum of an A0V star, HD 37887, observed near in 
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time and close in air mass to each object. 



2.2.1. Spectral Type Determination 



We list in Table [3] the spectral type of the host star of each of ou r transitional disks 



Most of our spectral- type information is gathered from the literature (Hillenbrand Il997b 



Rebull et al.l 12000 : lAllenl 1199a iDa Rio et al.1 120101 Fang et all 120091 : IParihar et al.l 120091 ) . or 
from unpublished results kindly provided by Lori Allen, John Tobin and Jesus Hernandez. 
In addition, we determined several new or improved spectral types from our SpeX spectra. 
In Figure [Tj we show the SpeX spectra of these objects and illustrate our spectral-typing 
procedure. We used Mg I, Al I, Na I, Ca I, and CO features in the H and K band as the 
SpeX sensitivity is best in those bands and as those absorption features are very sensitive to 
spectral types in the G-M range. 

It would have been best if the spectral resolution of our target spectra was the same 
as that of the standard spectra in order to distinguish adjacent lines and spectral depth. 
However, most of our SpeX spectra were taken at R ~ 800-1200. Only OriA-34 was observed 
with 0.3 arcsec slit giving R ~ 2000. Even though some lines are blended due to the modest 
spectral resolution, we were able to narrow down spectral types to about +/— 2-3 subtypes. 
We adopt the spectral type of the first five objects from top of Figure [T] from the spectral 
typing using our SpeX spectra. The S/N on OriA-302 is not good enough to determine a 
sub-class of its spectral type, but OriA-302 is thought to be an M-type star based on the 
broad spectral features attributable to Mg I, Na I, and CO overtone bands. We found that 
the uncertain spectral type from HECTOSPEC data (Lori Allen, private communication) of 
OriA-230, OriA-271, and OriA-294 are in reasonably good agreement with those from our 
spectra. 



2.3. Ancillary Data: Photometry 



We also compiled broadband photometry from ultraviolet to mid-infrared wavelengths 
from the literature. From the Naval Observatory Merged Astrometric Dataset (NOMAD) 
catalog ( jZacharias et al.ll2005l ) we collect B (0.44 /mi), V(0.55 / an), and R (0. 6 4 /mi) . There 



is a rather recent photometry data set for ONC objects from IDa Rio et al.l (120101 ) . From 



their tables, we gathered optical photometry at U (0.347 /mi), B (0.454 /im), V (0.538 /mi), 
TiO (0.6217 /mi), and I (0.862 /im). Photometry from one more optical band (I at 0.8 /im) 
was taken from the DENIS database as well as two near-IR bands, J (1.25 /im) and K (2.16 
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/xm). Most of our targets (except OriA-302) have 2MASS photometry in the J (1.25 /xm), H 



(1.65 /xm), and K (2.17 /xm) bands. We collected 2MASS information from ISkrutskie et al. 
( 120061 ). The photometry in the JHK bands for OriA-302 in Figure 2 are the averaged fluxes 
from our SpeX spectra. All of our targets were observed with IRAC (3.6, 4.5 , 5.8, and 8.0 
jttm) and MIPS (24 /xm) prior to the IRS observations (IMegeath et al.ll2012l ). The SEDs 
(Spectral Energy Distributions) constructed with the broadband photometry are shown in 
Figure 2. 



3. STELLAR PROPERTIES AND ACCRETION PROPERTIES 



3.1. Extinction Correction 



Extinction toward protoplanetary disks can lead to misclassification of evolutionary 
stages and misinterpretation of disk spectra. To minimize the effects of extinction toward 
our targets, we de-reddened our data based on the estimates of visual extinction (Ay) 
obtained using the following relationship between Ay and the color excess E(Xi — X 2 ) = 
([Ai] - [X 2 }) obs - ([Ax] - [A 2 ]) mt : 



Ai 



Ay 
Ax 2 



Ax 
Ax 



x E(X 1 -X 2 ) 



To get a color from two wavelengths, Ai and A 2 , we used 2MASS JHK photometry for most 
of the sample and DENIS IJH photometry for the rest, when available. We measured Ay in 
several ways for each object. We use either (J — H) 2 mass, (H — K) 2 mass, or (J — J)denis 
as an obser ved color ([Ai] — [A 2 1 ) ? ft.;. W e adopt I — J, J — H, and H — K photospheric 
colors from iKenyon fc Hartmannl (I1995T) o r J — H and H — K of the Classical T Tauri Star 
(CTTS) locus of colors from lMeyer et al.l (119971 ) as the intrinsic color, ([Ai] — [A 2 ]) int . To get 
A\, for each set of Ai and A 2 , we used the Mathis (1990) extinction curve for Ry = 5.0 if 
the resul t ing A y < 3. In case of Ay > 3, we followed the empirical extinction curves from 
McClurd (120091 ): two composite extinction curves, one for 3 < Ay < 8 and one for Ay > 8. 
After calculation of extinction corrections with each intrinsic color choice, we examined the 
extinction-corrected SEDs and selected a final result based on freedom from artifacts of the 
correction (e.g. artificial C0 2 ice features or structure in the silicate features) and good 
agreement with the photospheric spectrum of the star's type, at short wavelengths. (< 1 
/xm). In Table El we list the Ay method selected for each object: I — J, J — H, and H — K. 
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3.2. Stellar Properties 



We list the adopted spectral types in Table [3j Among 62 Orion A TDs, we have 55 
objects with well-determined spectral types, one with spectral type constrained to a broad 
range, and 6 with unknown spectral types. In Figure [3] we show the spectral type distribution 
of objects with known spectral types in ONC and L1641. The spectral types of ONC objects 
range from G5 to M5, while that of the L1641 objects are more concentrated around the Ml 
type. The SpT distribution of 27 TDs in ONC and 28 TDs in LI 641 generally ag ree with the 
SpT distribution of the general stellar populations of ONC (IRebull et al.l 120001 ) and L1641 
(IHsu et al.ll2012l ). respectively. 



The effective temperatures, T e jf, are adopted from iKenyon fc Hartmannl (119951 ) (see 
Table H|), corresponding to the spectral type of each object. For the objects of unknown 
spectral types, we used 3850K, the mean T e ff of Class II sources with known SpT in Orion A. 
The stellar luminosity (L+) of each object was derived from the stellar effective temperature 
and the stellar radius (R±). R+ was calculated from the scaling factor ((R*/d) 2 ) applied 
to the photosphere model, where d = 414 pc was assumed to be the distance to Orion A. 
The photosphere was derived from the intrinsic colors from lKenyon fc Hartmannl (119951 ) at 
temperature T e ff, scaled to match the de- reddened the 2MASS J b and photometry. The mass 
of star (M*) was inferred from the Siess PMS evolutionary tracks (ISiess et al.ll2000l ) using the 
assumed luminosities and effective temperatures. TDs with known spectral types in Orion A 
are shown on an H-R diagram along with Z=0.02 evolutionary tracks in FigureHl We also plot 
H-R diagrams for our TD sample from the Tau, Chal, Oph, and N1333 associations, which 
we also use in the present analysis, in Figure HJ We see that the host stars of transitional disks 
in Orion A lie furthest above the main sequence. This agrees with the generally-accepted 
age sequence in which Orion A is younger than that of Tau, Chal, or Oph. 

We compiled X-ray observations for our TD sample from a variety of sources. We 



searc hed for X-ray d ata in HE AS ARC, in the published literature (jGiidel et a. 



(2007) for 



Tau; Winston et al.l (120101 ) for NGC 1333), and in the Chandra Source Catalog (jEvans et al. 
20101 ) for objects in ONC. We wer e permitted pre-publi cation access to data from the XMM- 
Newton survey of L1641 (SOXS, iPillitteri et al.l 12013). We have also use d data from the 
Second XMM-Newton Serendipitous Source Catalog (jWatson et al.l l2009bl ) . The X-ray lu- 
minosity Lx we adopt in this work is that within 0.2-12 keV, the total band of XMM-Newton. 



These properties, T e jf, L+, M*, R+, and L x are listed in Table HI 
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3.3. Mass Accretion Rates: M 

We observed 52 of Orion A transitional disks with SpeX/IRTF in SXD mode from 0.8- 
2.4 jLtm to measure mass accretion rates from their hydrogen recombination lines: Pa7 (1.094 
jum), Pa/3 (1.282 /im), and Br7 (2.166 /zm). From the de-reddened SpeX spectra with the 
Ay determined as described above, we obtained mass accretion rates of all the objects except 
the G5 star, OriA-88, which shows strong Hydrogen absorption lines. 

The method of mass accretion rate measurement is as follows. We fit each hydrogen 
recombination line with a gaussian function plus a local continuum. We measure the line 
luminosity of each line from the fit. T hen we use the following relations to convert line 



luminosity to accretion luminosity L acc (IMuzerolle et al.lll998l ; iGatti et al.l 120081 ): 



log{L acc /L & ) = 1.36 x log(L Pa7 /L Q ) + 4.1 (2) 

log(L acc /L Q ) = 1.14 x log(L Pa p/L Q ) + 3.15 (3) 
log{L acc /L Q ) = 1.26 x log(L Bn /L Q ) + 4.43 (4) 



whence we obtain the disk-star accretion rate: 

In Figure [5] we show as an example the results for OriA-59. 

In general, the three recombination lines in our spectra yield similar results for accretion 
rate within a factor of 2-3, so we report the average in Table HI When fewer than three lines 
were detected we adopt the resulting average M as an upper limit and indicate them as such 
in Table H 



DISK PROPERTIES 



4.1. Transitional Disks: Selection Criteria and Their Variety 



Several spectral indices derived from IRS spectra of Class II YSOs have been used 
as a first step to separate transitional disks from the radially-continuous bulk of the pop- 



ulation (Fur 



Arnold et al. 



an et al.l 120091 : IWatson et al.l l2009al : iManoi et al.l 120111 McClure et al.l 12010 



20121 ). which we also use to identify the transitional disks in Orion A. The 



continuum spectral indices are defined as 

log(A 2 F A2 )-log(A 1 F Al ) 



log(A 2 ) - log(Ai 



(6) 



13 



The wavelengths, Ai and Ag, are selected to avoid emission features, and thus to represent 
the spectral shape of the optically-thick disk continuum emission. This in turn reveals disk 
structure, both radial (central clearings and gaps) and vertical (degree of flaring). The 
equivalent width of the 10 fim. silicate emission feature (EW(lOyum)), 

• 13/zm p _ p 

■T \ " A, con 



EW (10 fan) 



dX 



(7) 



8/im 



F\,con 

is a measure of the amount of optically thin dust per unit area of optically thick disk@ Thus, 
a large EW(lOyum) is a sign of a large amount of optically thin dust and/or an indication of 
a reduction in disk continuum due to the absence of optically thick disk in the region where 
the 10 micron emission feature is formed. Here we use the continuum spectral indices along 
with EW (10 fim) to identify disks with central clearing or radial gap. The principles are 
that small values of Uk-g — down to the color of photospheres — and large values of n 13 _ 31 
indicate the spectral "transition" that signifies a gap with outer gap radius in the few- to 
few-tens-of AU range; and that large values of EW (10 fim) connote warm optically thin dust 
in gaps. 

To find out break points for TDs among the Class II YSOs distribution in these observed 
parameter spaces, we utilized all (~600) IRS spectra of Class II YSOs observed in Orion 
A, T au, Chal, Oph, an d N133 3. Using the prop erties of TDs already well identified in 



Tau (IFurlan et al. 



( lArnold et al. 



20111 ). Cha I (IManoi et aUEOllt ). Oph (IMcClure et all 120101 ) . and N1333 



2012) clS cL guide, we found the breaks occur to at 



• tik-6 < -2.1 (the lowest octile for the distribution of njc-e) 

• "-13-31 > 0.5 (the highest octile for the distribution of ^13-31) 

• EW(10fim) > 4.3 (the highest octile for the distribution of EW (10 fim)) 



5 In our analysis EW(10/im) does not depend much on temperature and composition. As is evident in the 
ubiquity of 10 fim excesses and silicate emission features, disks around young stars or brown dwarfs always 
have a distribution of temperatures which exceed that required for efficient excitation of 10 fim continuum 
(T > 300 K), and therefore are well above that required for the silicate features, which are dominated by 
emission from cooler material at larger radii. We also know from the shape of the silicate emission profiles 
that essentially all the dust grain s we see are optically thin (internally, that is), and composed of amorphous 
and crystalline silicates (see, e.g.. lSargent et al.1 (|2009l )). In this case EW(10/xm) does not depend upon mass 
fractions of amorphous and crystalline material. One way to see this is to note that a grain with a given 
number N of silicate monomers has a fixed number of oscillators, with possibly a small range of oscillation 
frequency, in any given vibrational mode. Absorption or emission integrated over that mode is, to first order, 
proportional to N. (This is essentially the Thomas-Kuhn sum rule.) Thus two grains with the same N but 
different mineral fractions have the same equivalent width in the same vibrational mode - such as that which 
produces the 10 fim silicate feature - though the one with larger mineral fraction will have a larger number 
of oscillators at a small set of fixed frequencies, and identifiable sub-structure to the silicate feature. 



-14- 



In Figure [6] we plot ni 3 _ 3 i vs. nx-e, and si vs. £W(10/£m) for objects in Orion A. We 



also added objects in Tau (IFurlan et all 1201 ll ) in the plots for comparison 



To identify TDs in Orion A, we first selected objects satisfying one of the above condi- 
tions as possible candidates. Second, we rejected objects with spectral types earlier than G 
because they have significantly higher masses (> 2 M & ) and their circumstellar disks evolve 
much faster by possibly different disk clearing mechanisms than the case of transitional disks 
around the low-mass T Tauri stars. However, we keep objects with unknown spectral type 
as all appear to lie in the same luminosity range as the low-mass T Tauri stars. Third, 
we examined the detailed shape of the SEDs, and selected objects with deficits of infrared 
excess in the 2- 8 micron range, com pared to the appropriate median IRS spectrum of Class 



II disks in Tau (IFurlan et al.ll201ll ). We used the median of Taurus K5-M2 for objects with 
the spectral type earlier than M3, and the median of Taurus M3-M5 for objects with the 
spectral type of M3 or later. We adopt the Taurus median spectra instead of Orion A median 
because transitional disks in Tau and the median SED of Class II disks in Tau have been 
well studied and used widely for comparison in studies of other star-forming regions. As 
shown in Figure [7J the previously characterized TDs in Tau, Cha I, Oph, and N1333 pass 
these filters clearly. 

In Table we indicate how the TDs in this paper were selected based on the selection 
criteria. Three objects, OriA-44, OriA-88, and OriA-172, are classified as TDs based upon 
examination of their SEDs which are indicative of gap/central hole, despite the fact that 
they do not pass any of the three criteria. The spectra of 62 transitional disks in Orion A 
appear in Figure 2. We use these 62 TDs in Orion A and 43 TDs from Tau (13), Cha I (11), 
Oph (10), and NGC 1333 (9) to unveil properties of transitional disks in the following. 



4-1.1. Subclassification of Transitional Disks: CTD, WTD, and PTD 

The mid- infrared spectrum of a CTD like DM Tau and CoKu Tau/4, having a few- 
AU to few-tens of AU central clearing, is distinctive compared to a radially-continuous 
Class II disk: it shows very little continuum excess over the photosphere from near-IR (1- 
2 fim) to wavelengths in the mid-infrared (around 8-13 /im), at which point the excess 
increases exponentially with increasing wavelength until it matches or exceeds the median 
spectrum. The other two types of TDs have excesses from near-IR to the shorter mid-infrared 
wavelengths (~5-8 fim) that are smaller than or similar to the median. In some, like LkCa 
15 an d UX Tau A, veiling in the ne ar IR spectra points to an underlying optically thick inner 



disk (jEspaillat et al.ll2007bl . 120081 ). and their distinctive spectrum corresponds to optically 



thick inner and outer disks separated by a gap. These are PTDs. Intermediate between 
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CTDs and PTDs are WTDs, in which the weaker near- and mid-infrared excess is best 
explained by an optic ally-thin inner disk separated by a gap from an optically-thick outer 



disk, as in GM Aur (ICalvet et al.l 120051 ; lEspaillat et al.ll2010f ). Therefore, the distinction 



between CTDs, WTDs, and PTDs is whether an (optically thick/thin) inner disk exists or 
not. 

It is useful to define the Inner Disk Excess Fraction (IDEF) to characterize the near- 
infrared and shorter-wavelength mid-infrared excess fraction relative to the K5-M2 median 
spectrum of Class II sources in Tau. Since our sample is complete in H band data from the 
2MASS catalog and we have IRS spectra starting at ~5.2 /mi, we interpolate H band to 6 
/im to acquire the assumed spectrum covering 1.65-6 /xm: 

/•6/im _ 

Jl. 65 Um r \obj r \,photosphere u>A 

IDEI = fuTrn (°) 

J I 65/iUL A,median r X,photosphere u/ > 

In the case of OriA-88 for which we do not have an SL2 spectrum, the flux at IRAC channel 
3 (5.8 /im) is used instead of IRS fluxes. 

By taking the already well studied Taurus TDs as references, we adopt a set of infrared- 
excess ranges to subclassify TDs: 



CTD: IDEF < 0.25 
WTD: 0.25 < IDEF < 0.5 
PTD: IDEF > 0.5 



The IDEF values for the TD subclassification criteria are derived by adopting the K5- 
M2 median spectrum. The M3-M5 median spectrum is fainter than the K 5-M2 median 



spect rum due to an effect of lower stellar luminosity and lower disk emission (IFur lan et al. 



201 ll ). Therefore the excess emission over photosphere at 2-8 /im is weaker. It may lead 



to misclassification to compare objects of M3 or later spectral type to the K5-M2 median 
spectrum: an IDEF value of an object with M3 or later spectral type derived from the K5-M2 
median spectrum is generally lower than that derived from the M3-M5 median spectrum. 
Furthermore, there are no WTDs and few PTDs with M3 or later spectral type among 
the already well studied TDs in Tau, Cha I, and Oph. Thus, for M3-M5 types, we rely 
on confirmation by other TD selection criteria. Considering the additional contributions 
such as the disk inclination and scattered light to the degeneracy of the interpretation of 
the excess emission in the near-infrared, we consider the subclassification of TDs based on 
IDEF as preliminary. We note the subclassification based on IDEF values and some cases 
of exceptions of using this criteria in Table |5j 
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With these criteria we obtain 34 CTDs, 15 WTDs and 13 PTDs in Orion A. Adding the 
four other nearby associations of Tau, Chal, Oph, and N1333 brings the totals to 47 CTDs, 
17 WTDs and 41 PTDs. Spectral indices of these objects are plotted in Figures E] and [TJ 

Most CTDs are placed below the lower octile of riK-e, which reflect their negligible 
short-wavelength infrared excess from the inner disk, while PTDs are distributed over the 
whole range greater than the lower octile of n^-a- The WTDs are mostly located between 
the distribution of CTDs and PTDs in %_ 6 . 

In the plot of ni3_3i vs. EW(10fim) in Figure[6]and Figured we also indicate the ranges 
occupied by radially-continuous disk models with a range of inclination angles, stellar masses, 
degrees of dust settling, and mass accretion rates (ID'Alessio et al.ll2006l ; IWatson et al.ll2009ai 
Espaillatll2009l ). Most of TDs in Orion A lie outside the model polygon. However, in contrast 
to the positions of TDs in other star-forming regions, some of CTDs, WTDs, and PTDs in 
Orion A are located inside the polygon, indicating that they could be modeled as radially- 
continuous disks, with respect only to these two parameters. With the exception of IRS-18 
(PTD) and IRS-154 (WTD) in ONC, the transitional disks inside of the polygon are located 
in the region for vertically well-mixed disks (the area toward the upper right side in the 
polygon). Few obje cts lie in the domain of well-mixed disks for Tau, Chal, Oph and N1333. 
Fur lan et al.l (120091 ) concluded from this that substantial disk st ructural evolu t ion, e specially 
settling of dust to the disk midplane, has occurred in 1-2 Myr. Arnold et al.l (120121 ) showed 
that the disks in N1333 are statistically indistinguishable from those in Tau, Cha I and Oph 
in this regard. 

In the ternary plot of Figure [8] we see the distribution of TDs are separated from the 
radially-continuous disks in the three dimensional parameter spaces of uk-q, n i3-3ij and 
EW (10 fim): especially CTDs and WTDs are nicely located in the different region from the 
region occupied by the radially-continuous disks. 

We plot rii3_3i and EW (10 fim) of disks in Orion A along IDEF in Figure [9], as well as 
TDs in Tau as the references of criteria along IDEF. We find that ni^^i of TDs decreases as 
IDEF increase, i.e., CTDs tend to have larger n 13 _ 31 than PTDs. We note that CTDs with 
low IDEF span through the ranges of EW (10 fim), whereas most PTDs have high values 
(>4) of EW(10fim). 



4-1.2. TD Selection Criteria in the Literature 



We have compared our selection criteria for TDs to other selection criteria used in the 
literature. We plot TDs as well as Class II sources in OriA on the color-color diagrams 
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used for the selection criteria bvlMerin et all ( 120101 ) (left panel), ICieza et all (120101 ) (middle 



panel), and iMuzerolle et all ( 120101 ) (right panel) in Figure ITUl These authors used Spitzer 



IRAC and MIPS broadband photometry, not IRS spectra, for their selection criteria. 



Merin et al.l (120101 ) identified and characterized disks with inner holes from the Spitzer 



c2d Legacy program. They used [3.6]-[8.0] vs. [8.0]-[24] color-color diagram to identify TDs 
that fell in two separate regions: 

Region A: 0.0 < [3.6]-[8.0] < 1.1 and 3.2 < [8.0]- [24] < 5.3 
Region B: 1.1 < [3.6]-[8.0] < 1.8 and 3.2 < [8.0]- [24] < 5.3 



Region A selects TDs with central clearings like CTDs in our sample. Region B corre- 
sponds to the disks with some excess flux in the IRAC bands like WTDs and PTDs in our 
sample. Their selection criteria of region A agrees with our criteria for CTDs in L1641, but 
some number of CTDs in Orion A fall in Region B or outside the Region A and B (the left 
panel of Figure [TUJ) • 



Cieza et al.l ( 120101 ) set selection criteria based on the [3.6]- [24] vs. [3.6]- [4.5] color-color 
diagram, in which TDs are located in the region where [3.6]-[4.5] < 0.25 and [3. 6] -[24] > 
2. We see that most of our CTDs fall in that region, but PTDs have [3.6]- [4.5] > 0.25 (the 
middle panel of Figure [10]) . 



Muzerolle et al.l (120101 ) used spectral slopes (i.e. spectral indices) at two different wave- 



length intervals, 3.6-5.8 /im and 8-24 ju m. The criteria for the classical transitional disks 
with central holes based on a red as-24 i n IMuzerolle et al.l (120101 ) agrees with our criteria for 
CTDs. However, similar to ICieza et al.l (120101 ). most of our WTDs and all of our PTDs will 
be missed by this criterion (the right panel of Figure ITTJT) . 

From the comparisons of selection criteria for TDs, we find that (1) CTDs can be 
commonly identified by a variety of criteria; (2) a criterion utilizing th e largest wavelengt h 
interval in the IRAC channels (between chl (3.6 //m) and ch4 (8 //m) in lMerin et al.l (120101 )) 
can select WTDs and PTDs, but a color (i.e. the spectral slope) ba sed on other broadban d 



channels, such as [3.6]-[4.5] in ICieza et al.l ( 120101 ) or 3.6-5.8 /zm in IMuzerolle et al.l (120101 ). 
cannot find PTDs. 



4.2. Radial Properties of Transitional Disks: R wa ii 



One of the most important properties of a transitional disk is how large the gap or hole 
is; that is, the radius at which the inner wall of the optically thick outer disk lies. Bearing in 
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mind the limitations of inference of structure from the SED, we derive the radius, R wa ii, from 
the shape of the spectrum at the transition from small to large infrared excess within the 
IRS spectrum. O ur procedure, illustrated in Figure [LTJ is the same as in our previous work 
(IKim et al.ll2009l ): we model the inner edge of the outer disk as an optically thick insulating 
wall, with dust temperature T, as follows. We first subtract a power-law fit to the IRS 
spectrum in the 5-8 /im region to remove the flux from the photosphere or an excess from 
the inner disk. To the residuals, we fit a model with two components tightly constrained at 
13-16 /im and 30-33 /im. One component is emission from optically-thin astronomical-silicate 
dust ( IDraine fc Ledll984j ) with 0.1 /im radius, to represent the inner disk and the bulk of the 
optically thin atmosphere of the outer disk. The other component is a single-temperature 
(T) blackbody that represents th e insulated inner edge of the optically thick outer disk. As 
we described in lKim et al.l ( 120091 ). we do not aim to fit the details of the spectrum perfectly, 
but merely to separate the optically-thick continuum — for which the SED shape is the 
signature of the wall — from the silicate emission features centered at 10 and 20 /im. With 
the temperature T, we calculate R wa ii using radiative equilibrium at the inner wall: 



R 



wall 



4vraT 4 e /i? 



AuoT A 



(9) 



where A* is the effective albedo at stellar wavelengths and em is the effective emissivity at 
mid-IR wavelengths. We adopt (1 — A*)/ em = 1 : a perfectly black wall. 

We compared the resulting R wa ii of several transitional disks in Tau (CoKu Tau/4; 
DM Tau; GM Aur; LkCa 15; UX TauA), Cha I (CR Cha; SZ Cha; Til; T25; T35; T56), 
and Oph (Rox 44, 16126-2235AB) with the s ame quantity determined from detailed self- 
consistent models (jEspaillat et al.l I201CH . 1201 ll ). We found the derived R wa ii from Equation 
(9) are within about 33% of those obtained from detailed models. We also compared our 
Rwaii estimates to the cavity radius (R ca v) derived from the Submil limeter Array (SMA ) 
observation for DM Tau, GM Aur, LkCa 15, UX Tau A, and ROX 44 J Andrews et al.lboilh . 
and these results lie within 34% of each other. Comparison of R wa ii from the detailed model 
and the R cav observed by SMA also gives a similar range of difference, about 34%. Therefore, 
we adopt 33 % as the uncertainty of R wa u estimated from Equation (9). The R wa u results 
are listed in Table HI 



5. FREQUENCY AND AGE DISTRIBUTION 



Most previou s studies of transitional d i sks (e .g. IKim et al.ll2009l ; iMuzerolle et al.l 12010c 
Merin et all |201Q| ; ICurrie fc Sicilia-Aguilarl 120111 ) and protoplanetary disks (IFurlan et al. 
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20091 ; iLuhman et al.l l2010t iManojl 120101 ) have estimated the fractions of transitional disks 
to infer when the transition occurs and how long the transition stage persists. We take 
the fraction of CTDs+WTDs to compare to t he MuzerohVs CTDs b ecause the qualitative 
definition of the classical transitional disks in iMuzerolle et al.l ( 120101 ) is most nearly equiv- 
alent to the definition of the combined set of CTDs and WTDs in o ur sample. We note 
that the quantitative selection criteria used by IMuzerolle et al.l ( 120101 ) missed some WTDs 
as discussed in §4.1.2 

The fraction of TDs in this work is defined as n(TD type)/n(disks) using the number 
of TDs and number of Class II sources identified from the IRS survey of each region (see, 
Table |6]). We exclude samples in Ophiuchus for the frequency and age trend search even 
though it is listed in Table M because the selected TDs are from several sub-regions of 
Ophiuchus with differing estimated ages, and the sample sizes are too small to separate by 
subregion with high statistical significance. 



The TD fraction is plotted as a function of age in Figure [121 We note that (1) the 
fraction of CTDs+WTDs can be high over the broad age ranges (1-10 Myr) and (2) the 
fraction of TDs varies from region to region at young ages (< 3 Myr), with some having a 
fraction of only a few % while others have fractions of >20 % even at age < 1 Myr. 



Also shown in Figure [12] are the fractions of transitional disks of each type from this 
work only: CTD (square), WTD (circle), PTD (star), and the total of them (cross). If we 
consider the fraction of TDs including all CTD, WTD, and PTD, the fractions of all TDs 
in each region are significantly high even at very young ages of less than 1 Myr: 17±6% for 
N1333, 25±4% for ONC, and 26±5% for L1641 versus the TD fractions for the 2-3 Myr old 
associations: 8±2% for Tau and 16±5% for Cha I. 



6. DISTRIBUTION AND TRENDS OF PROPERTIES 

6.1. Transitional Disk Types vs. Spectral Type, R wa ii, and M 

In this section, we examine how the distributions of host star spectral types, R wa ii, 
and M compare for transitional disks with different inner disk structures: CTD, WTD and 
PTD types. In this analysis we exclude the TDs for which we do not have reliable host-star 
spectral types. 

In Figure [T3J we show the spectral type distribution for all TDs and for each subtype. 
The median spectral type differs slightly among the TD subtypes: M2 for CTDs; Ml for 
WTDs; K7 for PTDs; M0 for WTD+PTD; Ml for TDs. The spectral type distributions of 
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CTDs and PTDs differ noticeably, even if we consider the general spectral type uncertainty 
of one to two subtypes. However, a Kolmogorov-Smirnov (K-S) test shows that the spectral 
type difference between CTDs and PTDs is of marginal statistical significant: for CTDs vs. 
PTDs, D = 0.31 and p = 0.037; where D is the maximum deviation between the cumulative 
distribution of two groups and p indicates the probability that there is no significant difference 
between the distributions. 

In Figure [HJ we show the frequency distribution of R wa ii for each TD subtype. For 
almost 90% of TDs R wa ii is less than 30 AU. This pattern is similar whether the inner disk 
exists or there is an empty inner cavity. However, we see that the median R wa u is smaller 
for CTDs than WTDs and PTDs. Most CTDs have R waU < 10 AU (about 52 %), whereas 
71 % of WTDs and PTDs with inner disks have R wa ii > 10 AU. A K-S test that compares 
the Rwaii distributions of CTDs and PTDs yields values similar to those for the comparison 
of the spectral type distributions: D = 0.29 and p = 0.05, which is not very significant. The 
p value decreases to 0.02 when comparing CTDs and WTDs+PTDs. The difference in the 
R W aii distributions of the TD subtypes could also be due to the impossibility of distinguishing 
PTDs with very small gaps from radially continuous disks based on the IR spectrum alone. 
Overall, Figure HH shows that the IRS spectra are most sensitive to disk holes with R wa u < 
30 AU. 

In Figure [15] we plot the M distribution for all TDs and broken down by subtype. TDs 
with central clearings (CTDs) and gapped disks (WTDs and PTDs) differ significantly in 
M. Both groups also have M substantially smaller than the typical M of radially-continuous 
disks. The median M of the gapped disks is lO _8 ' 25 M /?/r and that of CTDs is 10~ 8 ' 7 
M & /yr. This visible difference is confirmed by statistical analysis with the K-S tests. When 
we include the upper limits, the K-S test results D = 0.46 and p = 0.001. Even when we do 
not include the upper limits, the statistical significant difference between the two groups is 
still valid with D = 0.5 and p = 0.002. Thus we confirm that the Ms of CTDs are smaller 
on average than those of WTDs and PTDs. 

For fair comparison of the M of radially-continuous disks and TDs, we need to be sure 
that these two groups have similar ages and ste llar masses because M decreases with ag e 
( IHartmann et al.l Il998[ ) and increases with M Q (IMuzerolle et al.l 120031 ; ICalvet et al.l 12004 ) . 
Because systematic diffe rences of stellar masses and ages can arise from a choice of different 
evolutionary tracks (e.g. ISimon et al.ll2000t iNajita et al.ll2007l ). the best would be to have M 
estimated with the same method and assumptions used for M estimation of TDs. However, 
there is currently no such M survey of low-mass T Tauri stars. A comparison between M 
of Orion A TDs and other r adially-con t inuou s Class II disks measured in an homogeneous 
method will be discussed by Kim et al.l ( 120131 . in preparation). 
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There fore, for the next 
Taurus by iGullbring et al.l j 



jest c omparison, we use the M measurements of T Tau ri stars in 



the evolutionary tracks of ID'Antona &: Mazzitel 



1998) for the following reaso ns. IGullbring et al.l (Il998h adopted 
il (|1997l) . which tends to result in younger 



stellar ages and lower s t ellar m asses. ISimon et al.l (l2000f ) found both the lBaraffe et al.l (119981 ) 
mode l and ISiess et al.l (12000 ) mode l agree with dynamical masses in the 0. 7-1 Mq range, 
while ID'Antona fc Mazzitellil (119971 ) do not agree as precisely. Even though IWhite fc Ghez 
( 1200 ll ) measured M of T Tauri stars in Taurus as ad opting an evolutionary track com- 
bined iBaraffe et al.l (119981 ) and iPalla fc Stahlerl ( 119991 ). the samples are binary systems 
and M of binary system is no t comparable to M of TDs. Hence, we adopt the median 
M from IGullbring et al.l ( 119981 ) as the median M of radially-continuous disks in Taurus, 
M = lO _7 ' 8 M /j/r@ While bearing in mind the possible uncertainties (~30%) of stel- 
lar masses b etween two different ev olutionary tracks, we confirm that the ages and stellar 
masses of the IGullbring et al.l (119981 ) sample are similar to those of TDs studied in this work. 
Thus our results confirm that TDs in general have substant i ally smaller M than rad ially- 
continuous disks, in a ccord with the findings by iNajita et al.l (120071 ) , iKim et al.l (120091 ) , and 



Espaillat et al.l (120121 ). 



6.2. Trends among TD Properties 

Understanding the correlations between disk and stellar properties of TDs is an essential 
and important key to understand how protoplanetary disks evolve from radially-continuous 
optically thick disks to a final planetary system. We have gathered the largest sample of 
TDs which are identified by homogeneous criteria. From this large sample, we are able to 
find not only trends of TDs in general but also detailed trends of the three different types 
of TDs. 

Except for M, we estimate the stellar/disk properties of TDs in Tau, Chal, Oph and 
N1333 in the same ways as for Orion A TDs. The M of TDs of other star-forming regions are 
from literature or personal communication. Adopting M measured by different techniques 
from that for the TDs of OriA should not affect significantly these trends, based on the 
demonstration of the tight correlation between the luminosities of IR hydrogen recombination 
lines and L nrn measure d from the blue excess spectrometry and/or U-band photometry by 



Muzerolle et al.l ( 119981 ). Therefore, we expect that any discrepancies in luminosity between 



6 We excluded three TDs, GK Tau, GM Aur, and IP Tau, in the sample of IGullbring et al.l (|1998h to 
measure the median M of radially-continuous disks in Taurus. If we include those three TDs, the median 
Mdot decreases to M = 10 8 M Q /yr. 
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the two (IR and UV) data sets on average will be small (IMuzerolle et al.l 119981 ). These 
stellar/disk properties and other properties obtained from the literature are also listed in 
Table H 

Among trends from many possible combinations of properties, we present the trends of 
interesting pairs of properties showing somewhat different behaviors from non-transitional 
disks in Figure [16] through Figure [2TJ The correlation parameters are calculated by using 
Unmix _err. pro which was de veloped for the Bayesian approach to linear regression with 



errors in both X and Y, by iKellyi (120071 ). and the trends line in each plot can be read 
in the manner of log 10 (Y) = (a ± e_a) + (/3 ± e_/3) log 10 (X). We indicate the correlation 
parameters (corr) and probabilities (P) of pairs of properties in Table [7] for TDs without 
separation by subtype and in Table El for two subgroups of TDs separated by their radial disk 
structures, i.e., CTDs and WTPs+PTDs. In general, a 5% or lower P value is considered to 
be statistically significant. Sometimes P < 2% is considered as a conservative threshold of 
statistically significant. Therefore, we interpret a correlation to be a statistically significant 
with P < 2% and to be a marginally significant if P is 2-5%. 



6.2.1. Trends of TDs 



Our search for trends related to M utilized weighted linear regression to account properly 
for upper limits of M in our samples. The resulting trend for M-M* in Figure [TBI (a) is M 
oc Mfc L6±0 " 3 . To com pare this to the M-M* relations fr om the previ ous studies (M oc M* 1 ,95 
jCalvet et al.ll2004r ): M cc M* 20 ( IMuzerolle etaHbooah : M oc M* 2 " 1 (IMuzerolle et~aHl2005h 



we include in Figure [TBI (a) a plot of the result, log M ~ 2 log M+ — 7.5 (IMuzerolle et al. 



2003 



Telleschi et al.ll2007l ). The slope of the M-M+ relation of TD host stars is roughly consistent 
with previous studies among T Tauri stars, but the trend line for M is shifted downward by 
factor of about 10 with respect to the thick dashed line representing T Tauri stars in Taurus, 
consistent with the results discussed in Section 6.1. 

In Figure fTBT-fb). we plo t Lx as a function o f M*. Compared to the results in Taurus 
\ogLx = 1.69 log M*+ 30.33 (ITelleschi et al.ll2007l ). the regression of our TD data is logL^ = 
1.1 log M± + 30.1: slightly smaller slope than that of T Tauri stars. 

Figure [TBI- (c) shows M as a function of Lx compared to the Taurus results derived from 
Taurus's M-M* and Lx-M*, M-Lx relation for T Tauri stars logM = 1.2 log Lx — 43. The 
M-Lx trend line for our TD sample, logM = 0.2 log Lx — 14.6. This shows that TD's M is 
not related to Lx, unlike CTTSs. 



The most interesting relations are the correlation of R wa ii to stellar properties. In Fig- 
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ure [T7T- (a) . we find a very strong correlation between M* and R wa ii- For both cases of M and 
L x , we also see the increasing trends of M-R wa u in Figure [TTl-(b) and L x -R W aU in Figure [T71- 
(c). However, we should recall that there is a tight correlation of M ic -R wa u and a statistically 
significant correlation of M-M*. Therefore, the correlation showing in the plots of M-R wa u 
may merely reflect the combination of correlations between M^-Rwau and M-M*. To test this 
hypothesis, we looked for a trend between M and R wa u in more restricted mass bins (e.g., 
0.2-0.4 M Q , 0.4-0.7 M , and 0.7-2.3 M ) and did not find significant trends. Similarly, the 
increasing trend of Lx-R W aii may be the result of combination of the correlations between 
L x -M* and M^-R waU . 

For further examination of the difference between M-R wa u relationship of TDs and that 
of CTTSs under the strong dependence on M±, we derived an expected M(M+)-R wa u of 
CTTSs by combining the M-M* of CTTSs (log M = 2 log M* - 7.5) and a strong M*-R wall 
correlation of TDs (logM* = OJlog-R^ - 1.0): log M(M*) = 1.4 log -9.5 (FigureO 

(b) ). The expected L X (M*)-R wa u of CTTSs considering the M* dependence is also derived 
by combining L x -M* of CTTSs and M*-R waU : log L X (M±) = 1.2 log Ruaii ~ 28.6 (FigureO 

(c) ). We find that the trends of TDs in M-R wa u and L x -R wa u generally agree with the 
trends of CTTSs. 

Several other interesting correlations for TDs are listed in Table The very strong 
correlation of L* ~ M^ 1 6 which is similar to that (L* ~ M* 1 ' 5 ) found from T Tauri stars 



in Tau by iTelleschi et al.l (120071 ) supports that the basic stellar properties, and L+, of 
TDs are not very different from the stars hosting radially-continuous flared disks and/or 
homologously evolved disks. 



6.2.2. Trends of M and L x for Transitional Disk Subtypes 

In this section we explore characteristics of M and L x , which are sensitive to the 
different inner disk structures and to the different disk evolution mechanisms, by searching 
for detailed trends of sub-samples grouped by different TD types. The correlation parameters 
of the detailed trend analysis for the case of inner clearings (CTDs) and the case of radial 
gaps (WTDs+PTDs) are listed in Table [HJ 

It is clear that the trend between M and M* of CTDs is different from that of WTDs+PTDs 
from the left panels ((a) and (b)) of Figure [TSJ When the inner disk is essentially empty as for 
the CTDs, no correlation exists between M and M*. However, for the case of WTDs+PTDs, 
M and M* are significantly correlated each other, and the relation (M oc M* 1 ' 9 ) is very close 
to that of CTTS (M oc M* 2 ). Therefore, the increasing tendency shown in Figure dH] prob- 
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ably leads to a strong effect from the significant correlation of M and M+ for WTDs+PTDs 
which still have an inner disk in their disk. In contrast to M, we find no significant difference 
between the two subgroups of TDs for the Lx-M+ correlation (the right panels ((c) and (d)) 
of Figure US]). 

Similar to the left panels of Figure [18] in the left panels of Figure [19] we also see no 
correlation of M-R wa u for CTDs but a significant and strong correlation for WTDs+PTDs. 
To find the genuine behavior of TDs between M and R wa u which is different from M of 
CTTSs, we divided M of our sample by the M(M*) (the thick long-dashed lines). The 
right panels of Figure [19] show the residual relation between M and R wa u of CTDs ((c); 
anti-correlation) and WTDs+PTDs ((d); no correlation) after taking the dominant effect of 
out. This supports the ideas that (1) the mass accretion behaves similarly to the case of 
radially-continuous disks while an inner disk exists in a TD, and (2) the mass accretion rate 
decreases as the size of inner cavity increases. 

In the left panels of Figure [201 CTDs have a very strong and significant correlation 
between Lx and R wa u, while WTDs+PTDs have a weak relationship with large uncertainty 
in the trend of Lx and R wa u comparing to that for all TDs combined. Therefore, we infer 
that the general increasing tendency of Lx along R wa ii in Figure [17] is due to the dominant 
effect of the strong correlation for CTDs. After removing the underlying contribution from 
M* to Lx-Rwaii relation by dividing Lx-Rwaii by Lx(M*)-R wa u, we see no correlation for 
CTDs and insignificant anti-correlation with large uncertainty for WTDs+PTDs from the 
right panel of Figure [2U] 

From Figure EH we confirm that M and Lx are not correlated with each other regardless 
of subtypes of TDs not like the strong correlation between them in case of T Tauri stars 
(the thick dashed line in the (a) and (b) panels). The residual M after removing the strong 
M-Lx correlation of T Tauri stars in the right panels are anti-correlated to Lx for both 
CTDs and WTDs+PTDs. We note that this is a similar to the trend found for CTTSs by 



Telleschi et al.l ( 120071 ). iDrake et al.l (120091 ) also shows the similar result that the objects with 
higher Lx have lower M and find a strong anti-correlation between Lx normalized to M* 
and M with CTTS/WTTS in ONC. Therefore, the anti-correlation shown in the right panels 
of Figure [21] for both CTDs (c) and WTDs+PTDs (d) may be a common characteristic of 
protoplanetary disks, not a unique characteristic of TDs. 
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6.3. Summary of Significant Trends 

Here, we summarize the correlations and trends of TD properties we find to help our 
insights on the origin of TDs, which will be discussed through the next section. 

L + and M* are independent of the subtypes of transitional disks. The significantly strong 
correlation between and R wa ii is consistent regardless of the subtype of TDs as M* oc 

R-wall ■ 

M and Lx vary/evolve with time. We have seen that some trends differ by subtypes of 
TDs and some trends are comparable to those of the larger T Tauri star population overall. 
The residual (or normalized) properties after removing the underlying effect of which is 
a parameter correlated with M and Lx may show the effect of the transition process from 
a radially continuous disk to a transitional disk. 

Some properties have very different trends according to whether the inner disk exists 
(WTD/PTD) or not (CTD). 

(1) Trends related to M show a strong correlation with WTDs+PTDs, but no correla- 
tion with CTDs: 

• M vs. M+: WTDs+PTDs show a very similar correlation to that of T Tauri stars with M 
oc M* 1 ' 9 . In stark contrast, no correlation shows for CTDs. 

• M vs. R wa ii- WTDs+PTDs show a very tight correlation between the two properties, M 
oc Rwaii 1 ' 9 , but the two properties are not correlated for CTDs. 

(2) Trends of some properties are strongly correlated for CTDs, but not for WTDs+PTDs: 

• Lx vs. M±: The trend of CTDs [Lx oc M* L8 ) alone is very close to the trend of T Tauri 
stars (L x oc M* 1 ' 7 ), but the slope (j3) of the trends of WTDs+PTDs is smaller by about a 
factor of two (L x oc M*). 

• L x vs. R wa ih A significant correlation between L x and R wa u for CTDs, but the trend for 
WTDs+PTDs is very uncertain. 

• residual M vs. R wa ii' The residual trend of M and R wa ii is opposite to the trend before 
removing M* effect from M vs. R wa iV- CTDs show an anti-correlation, but WTDs+PTDs 
do not show any correlation. 

On the other hand, M vs. Lx show no correlations from any subtype of TDs, in contrast 
to the strong correlation for T Tauri stars with full disks. 
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MECHANISMS FOR THE ORIGINS OF TRANSITIONAL DISKS AND 

CONSTRAINTS 



In this section, we briefly review proposed mechanisms for disk dispersal and compare 
how our findings are consistent, or not, with these predictions from the mechanisms: dust 
coagulation and settling; photoevaporation; inside-out disk clearing by MRI; gravitational 
effects of one or more low-mass companions. 

Grain growth and settling. Dust grains should grow and settle to the disk midplane 
during the protoplanetary disk evolution process. As grains grow larger, the opacity of the 
grains becomes smaller, leading to weaker continuum emission in IR range. On this basis, it 
has been proposed that the flux deficit shown in SEDs of TDs may be due to the existence 
of an opacity hole caused by grain gro wth and settling in the inner di sk rather than a rea l 



material deficit in the inner disk (e.g. iDullemond &: Dominikl 120051 ). iTanaka et al.l (120051 ) 



suggested that the SEDs of TDs could be due to differing opacity as a function of radius; 
e.g., sma ller disk optical depth in the inn e r disk than in the outer disk at 10 fim. Some 
authors (j Gar audi 120071 ; iBrauer et al.l 120071 . 120081 ) have demonstrated the short time scale 
of grain growth and settling. However, considering grain growth/settling only as the main 
mechanism of TDs cannot explain the sharp transition represented by the sharp edge of 
the inner wall of the outer disk which is implied from the distinctive c haracterstic of TP's 



SEDs and is support ed by many resolved submillimeter images of TDs (IHughes et al.l 12009 



Andrews et al 



201ll ).In any case, grain coagulation and settling is a process that takes 



monotonically longer at increasing radius within the disk, so it cannot possibly explain the 
gaps in WTDs and PTDs. 

MRI. A mechanism important in th e large scale mixing an d turb ulence in a disk is 
the Magnetorotational Instability (MRI). IChiang fc Murray- Clayl (120071) showed how MRI 
accele rates mass accretion and leads to inside-out disk clearing. iPerez-Becker fc Chiang 
(120 111 ) considered the ionization, necessary for the MRI mechanism, by stellar FUV radiation 
and demonstrated that the surface layer accretion driven by this could reproduce th e trend 
of increasing accretion rate with increasing hole size seen in TDs. ISuzuki et al.l ( ]2010l ). using 
MHD simulations with X-rays as the ionization source, showed the disk winds driven by 
MRI leads to a decrease in surface density in the manner of inside-out dispersal irrespective 
of the existence of a deadzone. 



These models based on the inside-out disk clearing by MRI draining predict correlations 
between R wa u, M, and L X - (1) a positive correlation between M and R wa ii at a given M* 
because a larger ionized area leads to more mass accretion; (2) a positive correlation between 
M and Lx at a constant R wa u because more ionization of disk surface as exposed by stronger 
energy sources leads to more mass accretion. However, our findings on the trends of M, Lx, 
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and R W aii do not support these predicted correlations between R wa u, M, and Lx- First, 
a strong correlation of M with R wa ii is shown only for the case of TDs with radial gaps 
(WTDs+PTDs) not for the TDs with inner cavities (CTDs). Furthermore, that correlation 
disappears when differing stellar masses are accounted for. Second, we found no correlation 
between M and Lx for TDs regardless of their subtypes. To test for correlation between M 
and Lx at a given R wa u, we examine trends for several different R wa ii ranges in Figure [221 
We find no correlation between M and L x for any ranges of R wa ii- 

Photoevaporation. Photoev aporation has long been thought to be an important and 



maior disk dispersal mechanism ( 


Shu et al. 


1994 


Font et al. 


2004; 


Alexander et al 


12006a 


b 


. Hi 



High energy radiation from the central star 
ionizes and heats the disk surface, which can be unbound and leaves as a wind beyond a 
certain radius. A gap can be opened when the accretion rate is small enough that ma- 
terial from the outer disk beyond the photoevaporation radius cannot replenish the inner 
disk. The inner disk, decoupled from the outer disk, drains on to the central star while the 
outer disk is evaporated into space. These models encounter contradictions when only the 
EUV ionizing radiation is considered: (1) the mass accretion rates of TDs is higher than 
the photoevaporation rate (or mass loss rate), (2) the ages of many TDs are less than the 
dispersal time scale of photoevaporation, and (3) the existence of the large radial gaps of 
PTDs. Recently, the photoevaporation model has been modified to address those issues by 



taking much higher energy radiation sources (FUV and X-ray) as sources ( jErcolano et al. 



2009: Gorti fc Hollenbach 2009; Owen et al 



2010. 



2011aU bl. 120121 ): (1) much higher mass loss 



rate can be driven by X- ravs (lErcolano et al.ll2009l ); (2) 1-10 AU gap creation at relatively 
early epochs (3-4 Myr) (iGorti et al.ll2009h . Considering d isk d ispersal by photoevapora- 
tion through EUV, FUV, and X-ray radiation, lOwen et al.l ( 120121 ) concluded that the stars' 
intrinsic X-ray luminosity should have a decisive role in disk's life times and evolution. 

We compar e our data to the t heoretical predictions from the X-ray photoevaporation 
(XPE) model by lOwen et al.l (120121 ) . Both the TD data and the XPE model have tendencies 
of increasing R wa u and M as Lx and M increase. However, most TDs of all three types 
in our sample do not fall into the model domains, as indicated in Figure 123^1 and Figure EU 
In particular, the XPE model fails to explain the objects with either large R wa u or the large 
M. More than half of all TDs fall outside the region allowed by the XPE in Figure 



Therefore, we conclude that X-ray photoevaporation is not the dominant mechanism for 
creating transitional disks. 



lOwen et al.l (|2012f) note the difficulty of determining a correlation between any two of the three properties 
because Lx, Rwaii, and M+ are interdependent. Therefore, the predictions in the Figures are the results 
from numerical simulation. 
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Morishimal ( 120121 ) developed a gas disk model taking into account layered accretion 
driven by MRI and X-ray photoevaporative winds. With a central star of 1 M & and an 
initial disk mass of 0.1 M , they found their gas dispersal model can open a gap at large 
radii while the mass accretion rate is still similar to that of CTTS when dead zones are 
considered. They show that their model with a dead zone can reconstruct the distribution 
of the observed transitional disks with high M and large R wa ii, which is the range X-ray 
photoevaporation cannot reconstruct. 

The average dis k mass of Class II obj e cts m e asured by observation s of th e submm-mm 
continuum of disks ([Andrews fe Williams! (120051 ) ; lAndrews Williams! (120071 ) ) in Tau and 
Oph is about 5 Mj. This is much less than the minimum mass solar nebular (MMSN) and 
the requirement of disk mass to form giant gas planets or multiple planets which are observa- 
tionally confirmed. It seems likel y that submillim e ter co n tinuum observati o ns sys tematically 



underestimate disk masses (e.g. iD'Alessio et al.l ( 120011 ). lHartmann et al.l (120061) ). Alterna 



tiyely, systems could have more substantial disk mass at the Class 0/1 stage (IGreaves fc Rice 
2010l ). 



Gravitational effects of companions. If the disks start with larger masses, photoevapora- 
tive gas dispersal models may encounter another challenge to deplete disk material because 
more massive disks requires stronger energy sources to evaporate disk material. On the other 
hand, for a more massive disk, the more favorable mechanism to open and clear a gap in a 
protoplanetary disk is giant planet formation by gravitational instability. Recent theoretical 
results on planet formation by gravitational instability suggest the possibili ty of gas giant 



planet formation at eve n much closer distan ces to the central star (10 AU by llnutsuka et al. 



( 120101 ); R < 25 AU by iMeru fc Batd (120 101 )). Even giant planet formation via core accre- 
tion model, which requi res longer time scales, is possible in ~1 Myr with dis ks substantially 
larger than the MMSN JPodson-Robinson et allbood ; breaves k Ricelboioh . 



Giant planet formation in proto planetary disks has been proposed as the origin of 



gaps / holes of TDs and debris disks (e.g. iHolland et al.lll998l ; lJura fc Turnerlll998t IZuckerman fc Song 



2004 ). Several hydrodynamical simulations show the companion's pres ence can reproduce 
the sharp inner disk truncation and gap formation in a shor t time scale ( jOuillen et al.l 12004 ; 
Rice et al.l[io06l : IVarniere et al]l2~00~6l Izhu et alJboilL l2012h . 



Whether the companion is a planetary object or a sub-stellar companion, the observa- 
tional features can be explained by the dynamical effects of the companion on the central 
star. The distribution of M shown in Figure US] strongly supports the idea of gap opening 
and disk dispersal by planet /companion formation. The displacements of M of the gapped 
disks and the CTDs from the median M of CTTS is almost a factor of 10, which is con- 
sistent with the estimated decrements of the mass flow across a gap created by a low-mass 
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companion jLubow fc D'AngelolliooS ). As we discussed in iKim et al.l ( 120091 ). the strong 
correlation between L+ or M+ vs. R wa u is reminiscent of the observed dependence between 
binary separation and the system's stellar mass. While this similarity might be taken to 
suggest that the stellar mass vs. separation dependence is imprinted when stars form, this 
argument has an important caveat. The strong trend of M+-R wa u of TDs is possibly driven 
by strong correlations of L*-M* and L±-R wa u. More detailed studies are needed to investigate 
this possibility. 

Giant planet formation can also explain the trends shown in M / M(M*)-R wa u (Fig- 
ure [19]). Assuming (1) M^isk oc M* at a given age which implies (2) normalization by M* is 
the same as normalization by Mdisk, the mass accretion rates of WTDs+PTDs may not be 
strongly related to R wa u since material accreting from an inner disk may not be dependent 
on a planet formation location or R wa u and the mass accretion from an inner disk may be 
dominant than mass accretion through a gap from an outer disk. In the case of CTDs when 
the inner disk is depleted, however, the outer disk mass will decrease as the size of an in- 
ner cavity increases. That will result less material to accrete through the larger inner hole, 
therefore the anti-correlation between M/M(M+) and R wa ii of CTDs is naturally explained. 

The disk-clearing com panion can be a stellar companion such as the case of CoKu Tau 
4 (jlreland &: Kraud 120081 ) because a large fraction of stars form in multiple systems. If a 
stellar companions are the dominant origin of TDs, they should be of lower mass than the 
primary to explain the strong correlations in M -M+ of WTD+PTD. Sensitive s e arche s have 
been made for 10 AU-sc ale binaries in Ta u rus (IPott et al.l (120101 ). Kraus et al.l (1201 lh ) . but 
only one (CoKu Tau/4; llreland fc Kraud (120081 )) of the nine transitional disk systems so 
studied has a stellar companion been detected with binary separation similar to the gap 
radius. The limi t s on c ompanion luminosity for the other eight rule out stars. In addition, 
Kraus fc Ireland! (120121 ) have tentatively detected an infant giant planet in the gap of the 
LkCa 15 disk. Therefore, planet formation could possibly be chiefly responsible for the origin 
of TDs. 



Currently it is not clear that how different paths of disk evolution are followed, and 
which mechanisms dominate on a track of evolution given different initial conditions. We 
do not deny the contribution to the dispersal of disk material and disk evolution from other 
mechanisms such as grain growth, MRI action in the inner disk, and photoevaporation. In 
other types of evolved disks, such as the anemic, homologously depleted, and weak excess 
disks, which are not covered in this study, grain growth and photoevaporation could be dom- 
inant mechanisms. However, we find that giant planet formation is probably the dominant 
mechanism in gap formation for our sample of TDs, based on the trends reported here. 
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8. Summary and Conclusions 

We presented the SEDs of 62 new TDs identified from IRS spectra in the Orion A star- 
forming region and discussed selection criteria for TDs. Utilizing the TDs already identified 
in Taurus, Chamaeleon I, Ophiuchus, and NGC 1333, in a manner similar to TDs in Orion 
A, we explored statistically robust trends with the largest and most homogeneous set of TDs 
to date. 

We presented a set of TD selection criteria and a quantitative empirical method to 
classify three subtypes of TDs: CTD, WTD, and PTD. 

We found the TD fraction of Class II YSOs is very high (~20%) even at the youngest 
ages (< 1 Myr, Orion A and NGC 1333). This could indicate early disk evolution even in 



the Class 0/1 stages, most likely due to giant planet formation ( ILubow fc D'Angeloll2006l ) 



We have examined various mechanisms of disk clearing utilizing our observation of 
TDs. We confirm that disk clearing mechanisms including MRI action and photoevaporation 
generated by X-ray/FUV photons are consistent with some of the observations, but not all. 
In particular, there are several observational trends which cannot be explained by these 
mechanisms, especially those related to M and Lx'- 

(1) the observed mass accretion rate suppression of CTD's (10~ 8,7 M & /yr) and PTD's (10~ 8 ' 25 
M Q /yr) compared to the radially continuous disks (10~ 7 ' 8 M & /yr); 

(2) the lack (or negative) correlation between M and R wa ii after accounting for M*'s effects 

(3) no correlation of Lx and R wa ii after correcting for M*'s effects 

(4) no correlation of L x and M at constant R wa ii- 

However, these properties/trends are naturally explained by substellar companions, formed 
within the disks in the early stages of disk evolution. Infant jovian-mass p l anets would 
nicely explain the size and structure observed for the gaps (IQuillen et all (12004 ): lEdear et all 



fl2007f n. 



Based on our results from the largest sample of transitional disks to date, as summarized 
above, we conclude that giant planet formation plays the dominant role in opening gaps and 
creating transitional disks presented in this work. 
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A. Notes on Individual Objects 



OriA-5: It is a single line spectroscopic binary (SB1), the average radial velocity is 
18±4.8 km/s, and the maximum velocity difference between two components is 3.8±4.8 
km/ s (ITobin et al.ll2009l ). Its SED shows steeply decreasing fluxes after 20 /im and this may 
be the reflection of the effect of a second companion. It also shows prominent crystalline 
silicate features. 



OriA-8: This one satisfies the criteria on both uk-g an d ^13-31- There are no obviously 
resolved sources within 30 arcsec around the target, but the target image in the Digitized 
Sky Survey (DSS) looks elongated, and it is suspected to have contributions from two sources. 

OriA-18: This object lies inside of the radially continuous disk model region in ri 13 _ 31 - 
EW(10/jm) space; however, its n K _ s is much less than the lower octile of n K _ 6 , which is one 
of the selection criteria. The reason for low ni3_3i is due to the decreasing flux after 20 /im, 
which could be the effect of external strong radiation evaporating the outer disk. There are 
3-4 other sources near this target. This object's spectral type is not known. 

OriA-39: It satisfies only one criterion of nx-6- Its IRS spectrum is very noisy (Fig- 
ure [2]). It is faint at IRS wavelengths, and the sky emission around it is very complicated 
because it lies in a fringe area of the bright nebula NGC 1977, even though there are no 
point sources within 40 arcsec. 



OriA-47: It lies in a very crowded and complicated region. There are many HH ob- 
jects about 2 arcmin away, and there is much complicated background emission from bright 
sources in the center of ONC, the Trapezium. There is also an additional point source about 
1.5 arcsec away. Recently, OriA-47 has also b een identified as a v ariable star, [PMD2009] 
185, as well as the nearby star [PMD2009] 183 JParihar et al.ll2009h . The IRS didn't resolve 
the signals from the two sources, and the IRS fluxes are the composite fluxes from both. 



OriA-88: Its SL2 data is not available because the array of SL2 is saturated due to the 
bright radiation entering in the IRS Peak-up cameras. This target is located just outside of 
the HII region, about 9 arcmin away from the center of ONC (M42). The PAH features are 
real and come from the disk. 
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OriA-149: This target lies in a complicated r egion of OMC2/3. It is reported to be an 
X-ray source ([TKK] 780) (ITsujimoto et al.l 120031 ). and there is a point source identified as 
an IR-source ([TKK] 774) about 5 arcsec away from the target. IRS SL and LL slits were 
placed to avoid [TKK] 774 as much as possible, but the source's IR radiation seems to affect 
LL1: there is a kind of extended emission entering in LL1 as background, very close to our 
target source. We tried to remove the effect from the extended emission as much as possible 
by blocking the pixels corresponding the emission and using multiple source extraction in 
AdOpt. 



OriA-154' Its spectral indices and EW (10 (Am) are similar to OriA-18, and it is lo- 
cated in a similar position on ni 3 _ 3 i-EW(10(Am) space. Source extraction for LL spectra 
was performed by AdOpt multiple source extraction because an additional source appeared 
in the LL slit. Its outer disk could be affected by strong radiation from nearby bright sources. 



OriA-164 '■ It lies in a dense dark core region, and it is very faint at optical wavelengths. 
This supports the large value of Ay found here, even though its spectral type is unknown 
so that there are large uncertainties in our Ay estimate. It is barely passed the criteria of 
EW(10(im). Its SED, which is restricted from J band to 35 (im in IRS data, with strong 
silicate features at 10 (tm and 20 (im, lets us infer that this object is a possible PTD. 



OriA-174'- The prominent emission lines in its IRS spectra are molecular hydrogen 
v = — > S(l), S(2), S(3) and S(5), which arise in the foreground and background cloud 
material rather than our target, and did not subtract away precisely. The object is an em- 
bedded source which is not shown up at optical wavelength ranges. 



OriA-188: The spectral type of this object is not known. There were no issues in the 
IRS data reduction, and its SED is like an object with a central clearing or gap. About in 
30 arcsec radius from the target there are three faint 2MASS objects, but they don't enter 
in the IRS slits. A significant environmental concern is that there are two OB type stars 
about 6.7 arcmin away from the target: iot Ori (09) and V2451 Ori (B7). The disk surface 
of the target may be affected by the strong radiation from the bright sources. 



OriA-198: This target also barely passed the criteria on EW (10 (Am) while its Uk-% and 
rii3_3i fail to pass the criteria. There is a bright nebula, [B77] 122, about 7 arcmin away, 
and no nearby sources within 10 arcsec of the target. OriA-198 has strong silicate features 



at 10 jum and 20 /im. These could be due to radial gaps in a more settled disk than a typical 
ClassII disk. 



OriA-221 Its EW (10) /im satisfies the TD selection criteria, and it lies in the area of the 
well studied TDs of Tau, the region of outliers in ni 3 _ 3 i vs. EW(10)fim space. After extinc- 
tion correction, this object's flux at optical wavelengths exceeds considerably that expected 
for photospheric spectra that fit well at longer wavelengths. This could be overcorrection of 
extinction toward the star, but we consider it more likely to be due to near-edge-on view 
and the extinction overcorrection of scattered visible light. There are several 2MASS objects 
within 40 arcsec, but they do not affect the source extraction at all. There is a bright B type 
star 7.4 arcmin away, HR 1911 (Bl; double or multiple system). This target is also a PTD 
candidate. 



B. Consideration of M and M w from XPE 

To ponder the mass dissipation throug h M and/or the ex pected mass loss rate (Mw), we 



overlay (Mw) by X-ray photoevaporation (jOwen et al.ll2012l ) on the observed mass accretion 



rate of TDs in Figure [23 The plus and cross data points in Figure [251 is the estimated Mw of 
each TD with its M* and Lx as the input condition of Mw- The estimated Mw have higher 
values than the mass accretion rates for most targets with high Lx and high M. However, 
Mw is weaker comparing M of CTTS indicated as the gray dashed line. From this, we may 
raise the question of when Mw overcomes M to create a gap/hole like in TDs. The M of 
TDs is not the M when Mw is able to open a gap. However, the estimated Mw are not 
very different whether for radially continuous disks (plus sign) or disks with a large inner 
hole (cross sign). If we assume that each TD should have followed the gray dashed line for 
their initial M-Lx relation before they became TDs, their initial M must have been greater 
than the theoretical Mw- And the current M of TDs are much lower than the theoretical 
Mw of TDs. Some possible scenarios might explain these inconsistencies: (1) the TDs in our 
sample have already evolved much after gap opening, and they have lost much inner material 
to accrete to the central star; (2) there could be some other contribution (e.g. planets) to 
make the mass accretion rate decrease quickly, such as intercepting disk material on the way 
to the central star while Mw forces material to drift away outward. 
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C. Trends of TDs in Orion A 

We examine trends for the sample of TDs in Orion A and the linear correlations of any 
two properties in log — log scale are listed in Table [9] and Table [10J Comparing trends in 
Table [9] and Table [7| before separating by subtypes of TDs, the trends of TDs in Orion A 
generally agree within 1 a uncertainties with the trends of the full sample of TDs not only 
from Orion A but also from Tau, Chal, Oph, and N1333, despite Lx and M* of TDs in 
Orion A tend to distribute toward lower values than that of other region. The trends of 
subtypes of TDs of Orion A in Table [10] are also not very different from that of full samples 
in Table M, even though uncertainties of the results from linear regression are larger due to 
smaller sample sizes for some correlations. We will discuss difference/similarity of properties 
of TDs and Class II obj ects in Orion A from that in other star-forming regions in a separate 



paper (IKim et al.ll2013l ) because the discussion is outside the scope of this paper. 
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Fig. 3. — Spectral type distribution for host stars of transitional disks in Orion A. 
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Fig. 4. — HR diagrams for host stars of transitional disks in this paper. Squares are for 
CTD s; circles are for WTDs; stars are for PTDs. Evolutionary tracks and isochrones are 
from lSiess et al.l ( 120001 ) (Z=0.02). Isochrone ages of various types of tran sitional-disk systems 
range from < 1 Myr to > 5 Myr. The average disk life time in Tau-Aur (IBertout et al.ll2007| ) 
is also shown dash-dotted line for reference. 



-55 - 




Fig. 5. — Example of mass accretion rate measurement from hydrogen recombination lines. 
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Fig. 8. — Disks samples of Orion A (black dots) and Tau (gray dots) are plotted in the spaces 
of three parameters. TDs are indicated by subtypes: CTDs (squares), WTDs (circles), and 
PTDs (triangles). Orion A TDs are open symbols and Tau TDs are filled symbols. The 
center of this plot is corresponding to the projected origin of the three axes. The three axes 
values are from to 100, which are linearly transformed from their original values of [-3, 
0] for uk-6, [-2, 3] for rii3_3i, and [-1, 9] for EW (10 fim). (A color version of this figure is 
available in the online journal.) 
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Fig. 11. — Example of fits using the simple model to estimate R wa ii of transitional disks 
in Orion A and other star-forming regions included in this paper. The thick solid line: 
a residuum of IRS spectra after subtracting a power-law fitted to 5-8 /im of IRS spectra, 
representing emission from the photosphere or a part of the inner disk. The dash-dotted 
line: astronomical silicate model to account for some contribution of emission from small 
dust grains in the atmosphere of wall and disk upper layers. The dashed line: a single 
blackbody profile with the wall temperature T. The dotted line: the continuum fit as the 
combination of the dash-dotted line and the dashed line. (A color version of this figure is 
available in the online journal.) 
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Fig. 12. — The fraction of transitional disks plotted against the estimated age of the associa- 
tion. The solid diamonds in the le ft panel are from the fraction of classical transitional disks 
defined by iMuzerolle et al.l (120101 ) . The empty diamonds in the left panel are the fraction of 
CTD plus WTD from this work. The right panel is for the shaded region covering ages of 
0-3 Myr in the left panel, and it shows fractions of TD (=CTDs+WTDs+PTDs) types in 
N1333, ONC, L1641, Tau, and Chal from younger ages to older ages of star-forming regions. 
The symbols in the right panel indicate the TD fraction (cross), the CTD fraction (square), 
the WTD fraction (circle), and the PTD fraction (star). 
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Fig. 13. — Spectral type distribution of transitional disks. The upper bigger panel shows the 
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Fig. 15. — M distributions. The upper panels include upper limits of M; the lower panels 
do not include upper limits of M. The left panels show M distributions separated by three 
different TD types; the middle panels compare the M distribution of disks with central 
clearings (CTD) and to those with gaps (WTD and PTD); the right panels result from 
adding all types of TDs. The results of K-S tests of M distributions between CTD and 
WTD+PTD in the middle panels: D = 0.46, p = 0.001 (upper middle); D = 0.5, p = 0.002 
(lower middle). 



- 66 - 



"D 
"D 



10"' 



10 s 



10 ! 



(a) 



♦ x 

+ ♦ 

I4f* * * * * 

. ♦ Xf± X X, 1 X 



io- 1 



IT 



10" 



0.1 



1.0 



m, (ivy 



(b) 



10' 



♦ 



10 3 



10 z 



x+ *x* 
x + ** iX 

«i* x + x 
F % x _♦ 



+ 



0.1 



1.0 



M. (Mg) 



"D 
"O 



10' 



10 "< 



10 s 



(c) 



♦ x» xX + 



>♦ + + 

+ 



X 

¥ 



10 



10" 



10 2 



10 30 
L x (erg/s) 



+ L1641 
XONC 
♦Tau 

Chal 
▼Oph 

N1333 

10 31 



Fig. 16. — Trends among M+, M and Lx- In each plot, different symbols indicate TDs in 
different star-forming regions (plus: L1641; cross: ONC; solid diamond: Tau; solid triangle: 
Chal; solid inverse-triangle: Oph; solid square: N1333). In each plot, the shaded area 
indicates the la uncertainty of the linear regression (of the logarithms). If a shaded area 
is narrow with high slope (e.g., M-M+), one can tell two properties in a panel is tightly 
correlated. If a shaded area is broad with very low slope (e.g., M-L x ), two properties in 
a panel is not correlated. The thick dashed li ne indicates a co r relatio n expected/observed 
among T Tau dis ks in Tau: M oc M* 2 from iMuzerolle et al.l ( 20031 ) in the upper panel; 
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Fig. 17. — Trends related to R wa ii with other stellar properties. Lines, shadow, and sym- 
bols have same meaning as in Figure [16j The thick long-dashed line indicates an expected 
trend derived from the relationship of a property in y-axis with M* shown in Figure [161 
and the strong M*-R wa u correlation shown in the panel (a): M oc Rwaii hi the panel (b); 
L x oc Rwaii 1 ' 2 the panel (c). (A color version of this figure is available in the online journal.) 
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Fig. 18. — Detailed Trends of M-M+ (the left panels) and L x -M+ (the right panels) separated 
by subtypes of TDs. The upper panels (a and c) show the correlation of CTDs (open square). 
The lower panels (b and d) show the correlation of WTDs+PTDs (solid star). The meanings 
of lines and shade are same as defined in the caption of Figure [161 (A color version of this 
figure is available in the online journal.) 
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Fig. 19. — Detailed trends of M-R wa u: correlations separated in two sub-groups, CTDs (open 
square) and WTDs+PTDs (solid star). The symbol in each panel is same as Figure [TBI The 
thick long-dashed line in the panel (a) and (b) is same in Figure [T71 (b). The right panels 
show the trend at no M* dependence by presenting the deviation of M from the thick long- 
dashed line of M(M+)-R wa u. (A color version of this figure is available in the online journal.) 
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Fig. 20. — Detailed trends of L x -R wa u- correlations separated in two sub-groups, CTDs 
(open square) and WTDs+PTDs (solid star). The symbol in each panel is same as in 
Figure fT8l The thick long-dashed line in the panel (a) and (b) is same as in Figure [T7] (c). 
The right panels show the trend at no M+ dependence by presenting the deviation of Lx 
from the thick long-dashed line of Lx(M±)-R waU . (A color version of this figure is available 
in the online journal.) 
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Fig. 21. — Detailed trends of M-Lx'- correlations separated in two sub-groups, CTDs (open 
square) and WTDs+PTDs (solid star). The symbol in each panel is same as in Figure fT8l 
The thick dashed line represent M-Lx relation of T Tauri star in Tau as explained in §6.2.1 
and Figure [TBI (A color version of this figure is available in the online journal.) 
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Fig. 22. — The tests to find the relationship between M and R wa u under a constant R wa u 
condition. We present the sub-groups separated by three R wa ii bins. The coverage of R wa u of 
each bin is indicated on each panel with the number of sub-sample (N), a linear correlation 
coefficient (corr) between logM and logLx, and a probability (p) of getting corr from random 
distribution. T^ enta u is Kendall's tau which indicates the degree of correlation between two 
variables; closer to 1, tighter correlation. p 2 indicate the two-sided p value of Tkentaii] if P2=l, 
the probability of no correlation is 100%. 
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Fig. 23. — Comparison of the predictions by X-ray photoevaporation (jOwen et al.l ( 1201 21 ). 
Fig 18 in their paper) to the properties of observed TDs. The gray shadows are the domains 
in which X-ray ph otoevaporation is dominant. The dashed line adopted from Fig 18. in 
Owen et al.l (120121 ) represents the maximum radius a TD may reach before thermal sweeping 
sets in and a disk dissipates beyond TD stage. The model domains ar e derived with R g oc 
and the positive correlation between and Lx ( jOwen et al.ll2012l ). similar to our finding, 
so the photoevaporation prediction and the observed TDs properties generally show positive 
correlations. 
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Fig. 24. — Co mparison of M vs. R m aii to the predictions by X-ray photoevaporation (gray 
shadow region: lOwen et al.l (120121 ) ) . Each symbol indicate different subtype of TDs: square 
(CTD); circle (WTD); star (PTD). 
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Fig. 25. — Comparison of M v s. of transitio nal disks to M w vs. Lx expected by the 
X-ray photoevaporation model (jOwen et al.ll2012l ). The plus signs are for the expected X-ray 
photoevaporation wind rates in case of a radially continuous primordial disk with TDs M* 
and Lx-, and the cross signs are for that in case of transitional disks with inner hole. The 
squares, circles, and stars are CTDs, WTDs, and PTDs of our sample. The gray dashed line 
indicate the expected M-Lx of CTTS, which is derived from M-M+ and L*-M* of CTTS in 
Figure [16j Detail discussion is in Appendix B. 
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Table 1 — Continued 



JN um. 


llxo name 


zivi/\oo name 


RA (J2000) 


DEC (J2000) 


A riRTrt 




Date observed 


region 


reduction 


W 


(2] 


(3] 




(4) 






(5) 




w 




\°) 


(9j 


(10) 


188 


8396035-5861 


05355047-0551422 


05 


35 


50 


.48 


-05 


51 


42 


3 


18838016 


44 


10/5/2007 


ONC 


AdOpt 


198 


8399269-6612 


05355825-0636431 


05 


35 


58 


.25 


-06 


36 


43. 


.2 


18798336 


39 


3/12/2007 


L1641 


man 


203 


8406280-6293 


05361506-0617369 


05 


36 


15 


.07 


-06 


17 


36 


.8 


18845696 


39 


3/12/2007 


L1641 


AdOpt 


204 


8406609-6247 


05361584-0614507 


05 


36 


15 


.86 


-06 


14 


50 


6 


18777088 


39 


3/12/2007 


L1641 


AdOpt 


211 


8411618-6426 


05362788-0625360 


05 


36 


27 


.88 


-06 


25 


35 


9 


18780160 


39 


3/12/2007 


L1641 


man 


218 


8416389-6503 


05363933-0630111 


05 


36 


39 


.33 


-06 


30 


11 


3 


18781184 


39 


3/12/2007 


L1641 


man 


219 


8416834-6225 


05364039-0613334 


05 


36 


40 


.40 


-06 


13 


33 


3 


18775808 


39 


3/9/2007 


L1641 


man 


221 


8416978-6185 


05364075-0611082 


05 


36 


40 


.75 


-06 


11 


08. 


.2 


18804992 


39 


3/11/2007 


L1641 


AdOpt 


223 


8424570-6484 


05365897-0629049 


05 


36 


58 


.97 


-06 


29 


04. 


.8 


18780416 


39 


3/12/2007 


L1641 


AdOpt 


227 


8440418-7404 


05373700-0724167 


05 


37 


37 


.00 


-07 


24 


16. 


6 


18814464 


39 


3/9/2007 


L1641 


AdOpt/auto d 


229 


8444779-6608 


05374746-0636298 


05 


37 


47 


.47 


-06 


36 


29. 


6 


18776064 


39 


3/12/2007 


L1641 


AdOpt 


230 


8445558-6860 


05374934-0651373 


05 


37 


49 


.34 


-06 


51 


37 


3 


18777600 


44 


10/6/2007 


L1641 


AdOpt 


237 


8457266-7161 


05381743-0709395 


05 


38 


17 


.44 


-07 


09 


39. 


.7 


18831360 


39 


3/9/2007 


L1641 


man 


239 


8464666-7838 


05383519-0750197 


05 


38 


35 


.20 


-07 


50 


19. 


.7 


18817536 


39 


3/9/2007 


L1641 


man 


257 


8502614-7795 


05400626-0747445 


05 


40 


06 


.27 


-07 


47 


44. 


6 


18779136 


44 


10/7/2007 


L1641 


man 


260 


8508501-7431 


05402040-0725540 


05 


40 


20 


.40 


-07 


25 


54, 


2 


18779392 


44 


10/7/2007 


L1641 


man 


270 


8519428-8120 


05404661-0807128 


05 


40 


46 


.63 


-08 


07 


12 


.8 


18821376 


36 


11/10/2006 


L1641 


man 


271 


8520557-7775 


05404931-0746327 


05 


40 


49 


.34 


-07 


46 


32 


6 


18777856 


44 


10/7/2007 


L1641 


man 


275 


8535781-7830 


05412586-0749506 


05 


41 


25 


.87 


-07 


49 


50 


9 


18806528 


36 


11/13/2006 


L1641 


man 


279 


8538911-7999 


05413338-0759562 


05 


41 


33 


.39 


-07 


59 


56 


5 


18820352 


36 


11/13/2006 


L1641 


man 


283 


8547523-7831 


05415405-0749534 


05 


41 


54 


.05 


-07 


49 


53 


.2 


18777344 


36 


11/13/2006 


L1641 


man 


290 


8562845-8151 


05423081-0809055 


05 


42 


30 


.83 


-08 


09 


05 


5 


18822912 


36 


11/13/2006 


L1641 


AdOpt 


291 


8564835-8250 


05423560-0815020 


05 


42 


35 


.60 


-08 


15 


01. 


.8 


18775296 


44 


10/7/2007 


L1641 


man 


294 


8567670-8803 


05424242-0848140 


05 


42 


42 


.41 


-08 


48 


13 


.8 


18830336 


36 


11/9/2006 


L1641 


man 


297 


8570949-8577 


05425027-0834378 


05 


42 


50 


.28 


-08 


34 


38. 


1 


18782208 


36 


11/9/2006 


L1641 


man 


301 


8576836-8303 


05430440-0818105 


05 


43 


04 


.41 


-08 


18 


11. 





18850304 


36 


11/15/2006 


L1641 


man 


302 


8579480-8565 


NA 


05 


43 


10 


.75 


-08 


33 


55 


.2 


18775552 


36 


11/9/2006 


L1641 


man 


303 


8580634-8516 


05431353-0831004 


05 


43 


13 


.52 


-08 


31 


00 





18776832 


39 


3/12/2007 


L1641 


man 



Note. — Column (1) Num.: the numbers are from the number sequence of 303 Class II objects observed in the IRS program 
number 30706. We use these numbers to identify objects easily. This number sequence will be used consistently in other 
future papers dealing with the objects in the PID 30706. 

Column (2) IRS name: sometimes IRS name indicated as in the column used to identify objects. 



Column (10) reduction: it indicate the methods of source extraction to get the SEDs in Figure 2. auto: an automated tapered 
column extraction in SMART with off-nod or off-order sky subtraction; man: a manual tapered column extraction in SMART 
with a polynomial sky subtraction; AdOpt: an optimal source extraction using an empirical point response function (PRF) 
in SMART; opse: an optimal source extraction using an analytical PRF. 

a SL & LL1: opsew7p; LL2: AdOpt. 

b SL: AdOpt; LL: opse. 

C SL: opse, LL:AdOpt. 

d SL: AdOpt; LL:auto onss. 
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Table 2. Observation log of SpeX spectra 



Num. 


IRS name 


obs. 


obs. date 


slit width 


total int. time comment 






semester 


(UT) 


(") 


(sec) 



1 


8336884-5290 


2011B 


11/6/2011 


0.5 


480 


5 


8343944-5609 


2011A 


2/28/2011 


0.8 


400 


11 


8347711-5533 


2010A 


2/18/2010 


0.8 


480 


13 


8350150-5595 


2010A 


2/17/2010 


0.8 


480 


16 


8353051-5229 


2010A 


2/18/2010 


0.8 


240 


19 


8355371-5480 


2010A 


2/17/2010 


0.8 


720 


23 


8357031-5072 


2011B 


11/7/2011 


0.8 


240 


24 


8358147-5505 


2011A 


2/27/2011 


0.8 


240 


25 


8358862-4842 


2010A 


2/16/2010 


0.5 


480 


26 


8359203-5026 


2011B 


11/10/2011 


0.8 


480 


29 


8361167-5475 


2010A 


2/18/2010 


0.8 


480 


34 


8365722-4816 


2010A 


2/16/2010 


0.3 


2520 


38 


8367284-5798 


2010A 


2/18/2010 


0.8 


240 


39 


8368137-4860 


2011B 


11/6/2011 


0.5 


480 


42 


8369106-5686 


2011B 


11/6/2011 


0.5 


720 


44 


8369980-5081 


2010A 


2/16/2010 


0.5 


240 


47 


8371984-5465 


2011B 


11/7/2011 


0.8 


480 


59 


8374388-6000 


2010A 


2/17/2010 


0.8 


480 


66 


8375158-5162 


2010A 


2/16/2010 


0.5 


960 


88 


8377167-5247 


2011A 


2/26/2011 


0.8 


120 


108 


8379312-5776 


2010A 


2/16/2010 


0.5 


240 


154 


8386523-5715 


2011A 


2/28/2011 


0.8 


320 


169 


8389033-4773 


2010A 


2/17/2010 


0.8 


480 


172 


8389673-4794 


2010A 


2/17/2010 


0.8 


480 


176 


8390285-5070 


2011B 


11/10/2011 


0.8 


240 


177 


8390554-6390 


2010A 


2/16/2010 


0.5 


240 


184 


8394216-5181 


2011B 


11/10/2011 


0.8 


400 


185 


8394454-5823 


2010A 


2/16/2010 


0.5 


240 


198 


8399269-6612 


2011A 


2/28/2011 


0.8 


240 


204 


8406609-6247 


2011A 


2/28/2011 


0.8 


960 


211 


8411618-6426 


2011B 


11/10/2011 


0.8 


240 


218 


8416389-6503 


2010A 


2/16/2010 


0.5 


480 


219 


8416834-6225 


2010A 


2/17/2010 


0.8 


480 


223 


8424570-6484 


2011A 


2/26/2011 


0.8 


960 


227 


8440418-7404 


2011A 


2/28/2011 


0.8 


360 


229 


8444779-6608 


2010A 


2/16/2010 


0.5 


240 


230 


8445558-6860 


2010A 


2/18/2010 


0.8 


720 


237 


8457266-7161 


2010A 


2/17/2010 


0.8 


480 


239 


8464666-7838 


2010A 


2/18/2010 


0.8 


240 


257 


8502614-7795 


2010A 


2/17/2010 


0.8 


480 


260 


8508501-7431 


2010A 


2/17/2010 


0.8 


240 


270 


8519428-8120 


2010A 


2/18/2010 


0.8 


240 


271 


8520557-7775 


2010A 


2/18/2010 


0.8 


960 


275 


8535781-7830 


2010A 


2/17/2010 


0.8 


240 


279 


8538911-7999 


2010A 


2/17/2010 


0.8 


960 
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Table 2 — Continued 



Num. 


IRS name 


obs. 


obs. date 


slit width 


total int. time comment 






semester 


CUT) 


(") 


(sec) 


283 


8547523-7831 


2010A 


2/17/2010 


0.8 


240 


290 


8562845-8151 


2011A 


2/28/2011 


0.8 


480 a, b, h 


291 


8564835-8250 


2010A 


2/17/2010 


0.8 


480 


294 


8567670-8803 


2010A 


2/18/2010 


0.8 


240 


301 


8576836-8303 


2011A 


2/28/2011 


0.8 


480 


302 


8579480-8565 


2010A 


2/18/2010 


0.8 


240 


303 


8580634-8516 


2010A 


2/18/2010 


0.8 


960 g 



Note. — (a) The source signal on order 8 is very faint, so that order 8 source 
extraction failed.; (b) It is a binary system with two sources resolved from SpeX 
observations. The possible primary source giving the adopted M in Table 0] is 
presented in the black line from 0.8-2.4 /im in Figure [5] The possible secondary 
is presented in the gray line in the same plot.; (c) The position angle of slit was 
27°.; (d) Resolved two sources on the guider image are dominantly displaced 
in the north-south direction with slight separation in east- west direction.; (e) 4 
among 6x120 sec exposures are combined in the spectrum.; (f) The two sources 
are aligned to N-S direction.; (g) The source signals in orders 6,7, and 8 is very 
faint so that source extraction for order 6-8 failed.; (h) The position angle of slit 
was -20°. 
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Table 3. Spectral Type and Extinction 



Num. 


IRS name 


Spectral type 


SpT reference 


A v 


Ay method 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


[ 


8*^6884 59Q0 


MO.O 


f 


0.61 


I- J 


5 


8343944-5609 

O O rfc O £7 'rl '-t ij \J \J iJ 


M5.5 


c f 


0.87 


l-.J 




83461 41-501 


M3.5 


f 


0.33 


I- J 


\l 


834771 1-5533 

OOrt I I -L J. OOOO 


K7.0 


SpcX 


8.27 


H-K 


13 


83501 50-5595 


M3.5 




o 


I- J 


16 


R^^ftM r >99Q 
oououu l-o^^u 


K1.0 


SpeX 


3.01 


J-H 


18 








\ 


CTTS 1 H 

\-J L J- O J 11 


19 


OOOOO f -L-O^iOU 


K7.0 


D 


0.15 


H-K 


23 


8^70^1-^079 

UOO 1 UO 1 uU i 


MO.O 


f 


1.88 


J-H 


24 


83581 47-5505 


K1.0 




4.21 


J-H 


25 


8358862-4842 


M3.5 


f 


1.62 


CTTS 1-H 


26 


8359203-5026 

UOOOiWO OUiiU 


M5.5 


f 


0.28 


I- J 


29 




MO.O 


SpeX 


0.73 


I- J 




8^(^799 481 fi 


Ml 


c nf ,Y 
opevv 


1 51 

1 .OL 


T T 
l-.l 


38 


8^fi798zL ^70,8 
oou f ^o^-o 1 yo 


MO.O 




1.1 


I-J 


39 


8^681 ^7-48fi0 


M3.5 


f 


0.65 


I- J 


42 


8369106-5686 


M4.0 


f 


1 


I-J 


44 


8^fiQQ80 ^081 
oouyyou ouo j. 


M1.0 


f 


0.71 


I-J 


47 


8^71 Q84. ^zLfi^i 


K7.5 


f 


4.68 


H-K 


53 


OO I 'J C/ 'J KJ U Ji'J 


M3.5 




0.24 


I-J 


59 


8374388-6000 


M2.5 


f 


0.28 


I-J 


66 


83751 58-51 62 


M5.0 


f 


0.22 


J-H 


88 


8377167-5247 


G5.0 


b 


2.57 


J-H 


108 


837931 2-5776 

OO 1 JOli O 1 1 \J 


KO.O 




1.82 


I-J 


149 


838536 5-51 06 






21.3 


CTTS 1-H 


154 


8386523-571 5 


M3.5 


e i' 

Si 1 


2.28 


H-K 


164 


8387753-4993 

OOO 1 1 iJJ 1C7i/J 






14.7 


CTTS 1-H 


169 


8389033-4773 


MO. 5 


f 


1.9 


H-K 


LIZ 


oooyO i o-4 ( y4 


t\/tq n 

IVlo.U 


I 




J-rl 


174 


8390073-5082 






8.42 


CTTS J-H 


176 


8390285-5070 


Ml. 5 


f 


0.79 


I-J 


177 


8390554-6390 


M2.5 


a 


0.12 


H-K 


184 


8394216-5181 


K7.5 


f 


1.79 


J-H 


185 


8394454-5823 


K6.0 


e 


1.57 


J-H 


188 


8396035-5861 






0.64 


CTTS J-H 


198 


8399269-6612 


K6.5 


a 


2.17 


J-H 


203 


8406280-6293 


K5.0 


d,i 


0.68 


CTTS J-H 


204 


8406609-6247 


M2.5 


d,f,h 


0.73 


I-J 


211 


8411618-6426 


M3.5 


a,f,h 


0.53 


I-J 


218 


8416389-6503 


M4.5 


a,f 


0.23 


I-J 


219 


8416834-6225 


M0.5 


e,f 


0.68 


I-J 


221 


8416978-6185 


K7.5 


a 


11.8 


H-K 


223 


8424570-6484 


M3.0 


e 


0.9 


I-J 


227 


8440418-7404 


MO.O 


a 


5.53 


CTTS J-H 


229 


8444779-6608 


MO.O 


e 


3.5 


H-K 
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Table 3 — Continued 



Num 


T t-i n m 




Qj_7 -L 1 LTI_>1 I^IIUL. 


4,. 


/lj/ Illl^ L11L/U. 






\ a ) 


w 


f51 
v'l 


y°> 


230 


8445558-6860 


MO.O 


a', SpeX 


2.42 


CTTS J-H 


237 


8457266-7161 


Ml. 5 


c 


4.59 


J-H 


239 


8464666-7838 


M1.0 


a 


0.39 


H-K 


257 


8502614-7795 


M3.0 


a 


0.63 


I-J 


260 


8508501-7431 


M2.0 


a 


1.17 


I-J 


270 


8519428-8120 


K7.0 


d,h 


9.53 


H-K 


271 


8520557-7775 


M2.0 


a', SpeX 


6.11 


H-K 


275 


8535781-7830 


M5.0 


c 


4.06 


H-K 


279 


8538911-7999 


Ml. 5 


h 


7.98 


H-K 


283 


8547523-7831 


M2.0 


c 


1.41 


J-H 


290 


8562845-8151 


M2.5 


a 


2.92 


l-.J 


291 


8564835-8250 


M2.5 


a 


2.81 


J-H 


294 


8567670-8803 


M2.5 


a', SpeX 


7.42 


H-K 


297 


8570949-8577 






19.9 


CTTS H-K 


301 


8576836-8303 


M3.5 


a 


3.16 


CTTS J-H 


302 


8579480-8565 


M 


SpeX 


3.41 


CTTS J-H 


303 


8580634-8516 


M3.0 


a,e 


6.56 


H-K 



Note. — Spectral type (column 3) and the methods or literature 
(column 4) from which we adopted the spectral type 

^Allen fc Mosbvl ||2008h . private communication; spectral types are 



measured from HECTOSPEC spectra 
a ' Allen fc Mosbv ( 20081) . private communication; spectral types 



measured from HECTOSPEC spectra, spectral types are highly un- 
certain 



feillenbrandl (jl997bh 
Rebull et all |200C| ) 



lAlleiJ (jl995l) 

Hernande j ( 2008 ) . private communication; spectral types measured 



from MDM spectroscopic data 



Hernandez fc Tobinl ([20091 ) . private communication; spectral types 



are measured from HECTOSPEC spectra 

f' 



Hernandez fc Tobinl (j2009l ): spectral types are measured from 



FAST spectra. 
tea Rio et all |2010l ) 
^Fang et al.l |2009l ) 
Parihar et all (|2009h 
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''it is M3.6 from iParihar et al.l (|2009l ). but for our analysis conve- 
nience we take it as M3.5 



Spex Spectral typing with SpeX spectra in this work 



Note. — Column (6) Ay method: the method which derived Ay in 
the column (5). See the text for details. 
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Table 4. Stellar and Disk properties 



Num 


IRS name 


T <=// 


U 




Rwall 


TD type 


M 


Lx 


region 






(K) 


(£©) 


(M e ) 


(AU) 




(M @ /yr) 


(erg/s) 




1 


8336884-5290 


3850 


0.64 





55 


13.7 


WTD 


(7.8±0.9)E-10 




ONC 


5 


8343944-5609 


3145 


0.54 





23 


3.0 


PTD 


<5.5E-10 




ONC 


8 


8346141-5010 


3420 


0.17 





31 


6.4 


CTD 






ONC 


11 


8347711-5533 


4060 


1.80 





70 


13.3 


CTD 


<2.37E-10 


(3.4±2.0)E+29 


ONC 


13 


8350150-5595 


3420 


0.45 





31 


20.0 


CTD 


(1.9±0.7)E-09 


(1.3±0.9)E+30 


ONC 


16 


8353051-5229 


5080 


3.65 


1 


88 


52.6 


CTD 


(1.6±0.1)E-09 




ONC 


18 


8355083-4835 


3850 


0.18 





58 


3.0 


PTD 






ONC 


19 


8355371-5480 


4060 


0.50 





77 


6.6 


CTD 


(3.2±0.6)E-09 


(2.1±1.0)E+30 


ONC 


23 


8357031-5072 


3850 


1.32 





54 


21.5 


WTD 


<1.38E-09 




ONC 


24 


8358147-5505 


5080 


5.98 


2 


19 


26.7 


PTD 


(3.6±1.2)E-08 


(2.1±0.5)E+30 


ONC 


25 


8358862-4842 


3420 


0.14 





30 


4.9 


WTD 


<2.3E-10 




ONC 


26 


8359203-5026 


3145 


0.17 





21 


6.6 


CTD 


(1.2±0.2)E-09 




ONC 


29 


8361167-5475 


3850 


0.41 





56 


10.3 


CTD 


(1.2±0.2)E-09 


(1.2±0.4)E+30 


ONC 


34 


8365722-4816 


3720 


0.25 





47 


7.6 


CTD 


<8.8E-11 


(2.0±2.5)E+29 


ONC 


38 


8367284-5798 


3850 


1.46 





53 


18.0 


CTD 


(1.7±0.6)E-09 


(4.0±0.7)E+30 


ONC 


39 


8368137-4860 


3420 


0.26 





31 


8.7 


CTD 


(7.2±4.8)E-10 




ONC 


42 


8369106-5686 


3370 


0.49 





30 


10.8 


WTD 


(2.6±2.4)E-09 


(4.7±3.0)E+29 


ONC 


44 


8369980-5081 


3720 


0.60 


1 


21 


5.9 


PTD 


(2.0±0.6)E-09 


(1.5±0.2)E+30 


ONC 


47 


8371984-5465 


3955 


4.37 





63 


50.6 


CTD 


<2.5E-09 


(2.6±0.4)E+30 


ONC 


53 


8373930-6325 


3420 


0.17 





31 


7.1 


CTD 






L1641 


59 


8374388-6000 


3525 


0.66 





36 


13.4 


CTD 


(6.7±0.9)E-09 


(1.8±0.2)E+30 


L1641 


66 


8375158-5162 


3240 


0.10 





23 


2.9 


WTD 


<2.3E-10 




ONC 


88 


8377167-5247 


5770 


62.34 


3 


10 


138.2 


CTD 




(1.0±0.1)E+31 


ONC 


108 


8379312-5776 


5250 


11.16 


2 


01 


59.8 


PTD 


(l.l±0.3)E-08 


(2.5±0.8)E+30 


ONC 


149 


8385365-5106 


3850 


7.40 





60 


61.2 


PTD 






ONC 


154 


8386523-5715 


3420 


0.86 





32 


14.1 


WTD 


(4.3±2.4)E-09 


(l.l±0.3)E+30 


ONC 


164 


8387753-4993 


3850 


5.02 





56 


18.5 


PTD 




(3.3±1.4)E+29 


ONC 


169 


8389033-4773 


3785 


0.81 





49 


7.9 


CTD 


(2.6±0.4)E-09 


(3.3±2.0)E+29 


ONC 


172 


8389673-4794 


3470 


0.62 





34 


7.2 


CTD 


<8.58E-10 


(3.7±0.3)E+30 


ONC 


174 


8390073-5082 


3850 


2.11 





53 


11.5 


CTD 






ONC 


176 


8390285-5070 


3650 


0.71 





42 


8.9 


WTD 


(3.5±0.7)E-09 


(7.9±1.2)E+29 


ONC 


177 


8390554-6390 


3525 


0.30 





36 


9.8 


CTD 


<l.lE-09 


(1.3±0.1)E+30 


L1641 


184 


8394216-5181 


3955 


1.12 





62 


16.7 


WTD 


(5.2±0.4)E-09 


(6.8±1.9)E+29 


ONC 


185 


8394454-5823 


4205 


1.30 





87 


11.7 


WTD 


(4.8±1.6)E-09 


(1.3±0.2)E+30 


ONC 


188 


8396035-5861 


3850 


1.19 





54 


9.3 


WTD 




(2.3±0.1)E+30 


ONC 


198 


8399269-6612 


4132 


1.33 





78 


19.4 


WTD 


(3.4±0.6)E-09 




L1641 


203 


8406280-6293 


4350 


1.46 


1 


06 


25.0 


PTD 




(3.3±0.1)E+30 


L1641 


204 


8406609-6247 


3525 


0.23 





36 


5.2 


CTD 


(1.5±0.8)E-09 


(8.3±4.2)E+28 


L1641 


211 


8411618-6426 


3420 


0.23 





31 


7.3 


CTD 


(5.7±4.7)E-10 


(2.0±0.5)E+29 


L1641 


218 


8416389-6503 


3305 


0.14 





26 


3.3 


CTD 


<4.73E-10 




L1641 


219 


8416834-6225 


3785 


0.50 





50 


11.4 


CTD 


(3.5±0.7)E-10 


(8.6±0.7)E+29 


L1641 


221 


8416978-6185 


3955 


1.12 





62 


17.8 


PTD 






L1641 


223 


8424570-6484 


3470 


0.34 





34 


9.6 


CTD 


(7.4±4.0)E-10 




L1641 


227 


8440418-7404 


3850 


0.81 





54 


11.3 


PTD 


(4.7±2.2)E-10 




L1641 


229 


8444779-6608 


3850 


0.61 





55 


13.7 


CTD 


(9.3±4.5)E-10 




L1641 
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Table 4 — Continued 



Num 


IRS name 


T eff 


£* 


Mi, 




TD type 


M 


Lx 


region 






(K) 


(£©) 


(M@) 


(AU) 




(M Q /yr) 


(erg/s) 




230 


844555o-oooU 


3850 


0.46 


0.56 


8.7 


CTD 


(l.y±U.2 JJi-lU 


(2.o±U.yjr J +zy 


L1641 


1 


04:0 / ZDD- / lul 


OOOU 


± .zo 


n aa 


1Q 1 
lO. 1 


O 1 lJ 


\\ ftTT" 1 n 

vO. IOJZj- 1U 




T 1 flzL1 

l^lO^il 


OQO 


Q/lft/lftftft 7CQO 

o4o4ooo- (ooo 


o ( zU 


1 Q1 

l.ol 


n a c 
U.40 


zO. < 


V V 1 JJ 


/i ftj_n ft Mr* ciq 
(i.oztu.o jJiz-Uo 




T 1 ft/1 1 
1j1041 


ZO / 


c^nofii /i 7"70c; 
oouzoi4- 1 ( yo 


q /i 7n 
o4 f U 


A OQ 

U.Zo 


n qq 
U.oo 


A 1 
4. 1 


L 1 L) 


/ 1 QJ_1 n\T7 no 

(i.ozti.u jji/-uy 


(z.u±z.ujJi/+zy 


T 1 ft/1 1 
1j1041 


OftO 
ZOU 


Qc;nctrni 7/iQi 
OOUoOUl- / 4ol 


0O0U 


n Qfi 
U.oo 


n qq 
U.oo 


n i 

y. i 


L 1 L) 


/q /ij_n ft Mr* no 
(o.4±U.O JH/-Uy 


(i.o±i.4jJi/+zy 


T 1 ft/1 1 
1j1041 


z /U 


o0iy4Zo-olzU 


/i nftn 
4U0U 


O 7Q 

y. to 


n qq 
U.oo 


oy.o 


L 1 L) 


c fti "c 1 in 
<0.OlJi-lU 


(a n_i_n c^^T7 , i Qn 
(4.U±U.0 JJI/+0U 


T 1 ft/1 1 
1j1041 


ZY 1 


OOZUOO (-1(10 


0O0U 


U.4 / 


n qq 
U.oo 


1 q n 

io.y 


L 1 iJ 


^ 1 ft/i t? in 
<C1.04J±/-1U 


/1 ftj_n n\T? 1 00 

(i.o±u.y jji/+zy 


T 1 ft/1 1 
L;l041 


ZYO 


oOoO <ol- tool) 


OZ4U 


U.oO 


n oc; 
U.zO 


y.o 


W 1 JJ 


^q ,1 t? no 
<o.o4Ji-Uy 


/o c;_i_n oMr 1 Qn 
(z.0±U.z JJ±/-|-oU 


T 1 ft/1 1 
1j1041 


97Q 


Rt^SQI 1 7QQQ 

ooooy i i- ( yyy 


OOOU 


i .uo 


n 41 


1 ^ Q 

±o.y 


VV 1 LJ 


^9 1 ^TT 1 n 

<^Z. IOJZj- 1U 




ijlO^tl 


zoo 


QK/IT'KOQ 7CQ1 

004 1 OZO- / Ool 


0O0U 


n e; q 
U.Oo 


n qq 
U.oo 


1 Q 1 




/ q 1 i n n\T? no 

(o.iztu.y jji/-uy 


/ft c;j_n o\t? 1 Qn 
(o.0±U.z JJI/+0U 


T 1 ft/1 1 
1j1041 


zyu 


fic:ft9fi/1 £ Q1 c;i 
000Z040-O101 


OOZO 


n /in 
U. J.U 


n Qft 
U.oo 


A 1 
4.1 


r 1 IJ 


9 Q-U9 a Mr no 
^z.o±z.4j Ji/-uy 




T 1 ft/1 1 
1j1041 


zy l 


oOo4ooO-ozOU 


OOZO 


n a q 
U.4o 


n Qft 
U.oo 


Q O 

o.z 


L 1 iJ 


/"o n_Lo q\t? no 
(z.UztZ.o JJl/-Uy 


/1 Qj_n ft Mr* 1 on 
(l.o±U.o JJi+zy 


T 1 ft/1 1 
1j1041 


zy4 


QKftTft'T'n OOflQ 

oOD / /U-ooUo 


Q c:oc; 
oOzo 


z.UU 


n or 
U.oO 


OQ Q 
ZO.O 


L 1 iJ 


/ 1 oot? no 
<Cl.zZJl/-Uy 


(a kxq ftMr 1 00 
(4.0±o.o JJl/+zy 


T 1 ft/1 1 
1j1041 


zy / 


Q^YOO/IO Q^'7'7 

oo / uy4y-oo I ( 


ooOU 


1 7^ 


n c;q 
U.Oo 


99 7 
ZZ. ( 


r 1 IJ 






T 1 ft/1 1 
1j1041 


QOI 

oUl 


Q^VftQQft QQm 

OO / OooO-ooUo 


q /i on 
o4ZU 


n qq 
U.oo 


n qi 
U.ol 


1/19 
14. Z 


I > T 1 "> 
r 1U 


z' Q-i-Q 1 Mr in 
^O.ozto.l JJi-lU 


i ( o.o±i.uji!j+zy 


T 1 ft/1 1 
1j1041 


Qno 
oUz 


OO / y4ou-oooo 


ooOU 


n on 
U.ZU 


n c; Q 
U.Oo 


0. f 


W 1 JJ 


(O.oztl. f JJ1/-1U 




T 1 ft/1 1 
L;l041 


QnQ 
oUo 


OOoUOo4-o010 


Q/1 7H 
o4 ( U 


n /i q 
U.4o 


n qq 
U.oo 


1 n /i 

1U.4 


^ 1 JJ 


/9 O-Un oW no 
i^z.yztu.o jji-uy 


( 1 7-un sWj-Oq 
^1. 1 ±u.o ji!j±zy 


T 1 ft/1 1 
1j1041 




U4zUzoUo 


3420 


0.19 


0.31 


12.3 


PTD 






Tau 




U41z5+zyUz 


3685 


0.36 


0.44 


17.0 


CTD 






Tau 




Ooivu lau4 


3650 


0.60 


0.44 


9.4 


O 1 JJ 


^- 1 |"iT7 1 na 
<l.UJi-lU 




Tau 




DM lau 


3720 


0.20 


0.44 


3.9 


U 1 JJ 


/q 1 j_i 1 Mr nob 
(o.l±l.l ;J1/-Uy 




Tau 




GKTau 


4060 


1.37 


0.71 


9.7 


PTD 


/ft n_i_o 1 \ T7 1 nnc 

(o.U±z.ljJi/-Uy 


/o 1 n \ tt 1 1 nn 

(8.8±0.2jhj+29 


Tau 






4350 


1.00 


1.25 


26.9 


W 1 JJ 


/ /i 7*_L1 ft\T7 nnb 

(4. 7 ±1.0 j-hj-uy 




Tau 




HKTau 


3785 


0.65 


0.50 


16.1 


PTD 


(2.z±0.ojJi/-Uo 


(0.1±Z.ljJl/±ZO 


Tau 




T "D rp 

IF lau 


3850 


0.69 


0.55 


12.6 


r 1 JJ 


1a oj_i a\t? 1 nb 
(4.Z±1.4JJ1/-1U 




lau 




LkCal5 


4350 


1.00 


1.08 


47. 


PTD 


(3.3±l.l)hj-09 




Tau 






a nftn 
4U0U 


n 07 
U.o I 


n *7Q 
U. / o 


ft n 

o.y 


I >T 1 1 

r 1 JJ 




/I ftj_n 1 ~\T7 1 OO 

(i.o±u.i jji/+zy 


Tau 




RYTau 


5945 


19.75 


2.36 


27.6 


PTD 


/ 1 I 1 \ tt 1 nob 

(9.1±3.1jhj-08 


/o n 1 n iXtt 1 i on 

(o.U±U.ljJi+oU 


Tau 




UaI auA 


4900 


2.20 


1.63 


59.8 


PTD 


/1 1 ^ \ T7 nob 
(l.l±0.4jtj-08 




Tau 




\T A 1 fl Wot rft 

v 4iu Arayo 


Q1 A X. 
O140 


n ack 

U.4D 


n 99 
u.zz 


7 Q 

i .y 


o i jj 




(1 Q-U1 fiW_LOQ 
( .O±1.0 )S2j-\-Zo 


Ton 

lau 




2MJ1124118 


3240 


0.06 


0.20 


2.9 


CTD 






Chal 




Baud43 


3445 


0.25 


0.32 


6.5 


PTD 




(5.6±0.3)E+28 


Chal 


T54 


CHX22 


5520 


4.10 


1.66 


37.1 


CTD 






Chal 




CHXR22E 


3420 


0.20 


0.31 


7.4 


CTD 




(2.6±0.3)E+29 


Chal 




CRCha 


4900 


2.54 


1.70 


16.9 


PTD 


(8.8±3.0)E-09 b 


(1.6±0.0)E+30 


Chal 


Til 


CSCha 


4205 


1.00 


0.87 


41.1 


CTD 


(1.2±0.4)E-08 b 


(2.8±0.1)E+30 


Chal 




IS091 


3470 


0.30 


0.33 


2.2 


CTD 




(2.8±1.0)E+28 


Chal 


T25 


Szl8 


3470 


0.20 


0.33 


8.1 


CTD 


(9.7±3.3)E-10 b 




Chal 


T35 


Sz27 


3850 


0.40 


0.56 


15.3 


PTD 


(1.2±0.4)E-09 b 


(5.8±5.8)E+29 


Chal 


T56 


Sz45 


3785 


0.30 


0.52 


18.3 


WTD 


(1.5±0.5)E-09 b 




Chal 


T6 


SZCha 


5250 


2.40 


1.55 


29.5 


PTD 


(2.4±0.1)E-09 b 


(1.2±0.1)E+30 


Chal 




16126-2235 


3470 


0.80 


0.34 


12.2 


PTD 




(1.9±0.0)E+30 


Oph off 




16201-2410 


6030 


4.90 


1.41 


29.2 


PTD 






Oph off 




16220-2452 


3470 


0.28 


0.33 


8.7 


PTD 






Oph core 




16225-2607 


4060 


1.37 


0.71 


14.2 


PTD 






Oph off 



Table 4 — Continued 



Num 


IRS name 


T e ff 








TD type 


M 


Lx 


region 






(K) 


(£©) 


(AT©) 


(AU) 




(M Q /yr) 


{erg/s) 






16289-2457 


5770 


4.61 


1.53 


24.7 


PTD 






Oph off 




DOAR28 


4350 


1.00 


1.08 


15.0 


CTD 






Oph off 




GY195 


3470 


0.38 


0.34 


6.6 


PTD 


(1.0±0.5)E-09 d 




Oph core 




IRS3 


3580 


0.39 


0.38 


13.4 


PTD 


(1.2±0.4)E-08 d 


(5.0±1.2)E+28 


Oph core 




Rox44 


4900 


3.00 


1.77 


26.7 


PTD 


(9.0±3.1)E-09 b 


(6.1±0.3)E+29 


Oph off 




SR9 


4350 


2.03 


1.06 


14.6 


PTD 


(1.2±0.4)E-08 




Oph core 


spID 50 


LAL 106 


3309 


1.63 


0.27 


45.4 


PTD 




(1.2±0.3)E+29 


N1333 


spID 52 


LAL 110 


4060 


1.01 


0.71 


18.0 


PTD 




(2.2±0.1)E+30 


N1333 


spID 57 


LAL 129 


4060 


0.92 


0.72 


23.2 


PTD 




(3.7±1.6)E+28 


N1333 


spID 73 


LAL 171 


4900 


8.38 


2.45 


32.2 


PTD 




(8.4±0.2)E+30 


N1333 


spID 101 


LAL 245 


3370 


0.07 


0.26 


8.4 


PTD 






N1333 


spID 110 


LAL 279 


3525 


0.02 


0.31 


3.7 


PTD 






N1333 


spID 116 


LAL 300 


4205 


1.17 


0.87 


5.0 


PTD 






N1333 


spID 136 


LAL 321 


3955 


0.83 


0.63 


12.3 


CTD 




(1.7±0.1)E+30 


N1333 


spID 137 


LAL 331 


3955 


0.23 


0.66 


6.6 


CTD 




(1.3±0.2)E+29 


N1333 



Note. — The calculation of R wa ii is explained in Section 4.2, and the values listed in the column have 
an uncertainty of about 33 %. M of ONC and L1641 is from this work. M for other regions are from 
the literatures. Lx are all from XMM- Newton Serendipitous Source Catalog (jWatson et al.ll2009bf ). 
The energy band is total energy band of XMM-Newton (0.2-12 keV). 



Naiita et all (|2007h 



Espaillat et al.l (|201ll ) 



Giidel et alJ (120071 ) 
( 2006h 



Natta et al 
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Table 5. Sub-classification of TDs 



LLJ 


name 


region 


n-K-6 


"13-31 


cj vv ^iu^/m ) 


subtype 






( 2 ) 


( 6 ) 




\?) 


W 


V) 


K°) 


1 


8336884-5290 


ONC 


N 


N 


Y 


WTD 


0.43 


5 


8343944-5609 


ONC 


N 


N 


Y 


PTD 


0.52 f 


8 


8346141-5010 


ONC 


Y 


Y 


N 


CTD 


0.27 c 


11 


8347711-5533 


ONC 


Y 


Y 


N 


CTD 


0.00 


13 


8350150-5595 


ONC 


Y 


Y 


N 


CTD 


0.06 


16 


8353051-5229 


ONC 


Y 


Y 


N 


CTD 


0.07 


18 


8355083-4835 


ONC 


Y 


N 


N 


PTD 


0.58 


19 


8355371-5480 


ONC 


Y 


Y 


N 


CTD 


0.05 


23 


8357031-5072 


ONC 


N 


N 


Y 


WTD 


0.51 d 


24 


8358147-5505 


ONC 


N 


N 


Y 


PTD 


0.62 


25 


8358862-4842 


ONC 


N 


Y 


Y 


WTD 


0.39 


26 


8359203-5026 


ONC 


Y 


N 


N 


CTD 


0.13 


29 


8361167-5475 


ONC 


Y 


Y 


N 


CTD 


0.11 


34 


8365722-4816 


ONC 


Y 


Y 


Y 


CTD 


0.18 


38 


8367284-5798 


ONC 


Y 


Y 


Y 


CTD 


0.07 


39 


8368137-4860 


ONC 


Y 


N 


N 


CTD 


0.18 


42 


8369106-5686 


ONC 


Y 


N 


Y 


WTD 


0.41 


44 


8369980-5081 


ONC 


N 


N 


N 


PTD a 


0.65 


47 


8371984-5465 


ONC 


Y 


Y 


N 


CTD 


-0.78 


53 


8373930-6325 


L1641 


Y 


N 


N 


CTD 


0.15 


59 


8374388-6000 


L1641 


Y 


Y 


N 


CTD 


0.22 


66 


8375158-5162 


ONC 


N 


Y 


N 


WTD 


0.40 


88 


8377167-5247 


ONC 


n/a 


n/a 


n/a 


CTD b 


0.22 


108 


8379312-5776 


ONC 


N 


N 


Y 


PTD 


0.90 


149 


8385365-5106 


ONC 


N 


Y 


N 


PTD 


0.12 cl 


154 


8386523-5715 


ONC 


Y 


N 


N 


WTD 


0.33 


164 


8387753-4993 


ONC 


N 


N 


Y 


PTD 


0.28 cl 


169 


8389033-4773 


ONC 


Y 


N 


N 


CTD 


0.13 


172 


8389673-4794 


ONC 


N 


N 


N 


CTD 


0.36 c 


174 


8390073-5082 


ONC 


Y 


N 


N 


CTD 


-0.33 


176 


8390285-5070 


ONC 


N 


N 


N 


WTD a 


0.41 


177 


8390554-6390 


L1641 


Y 


Y 


N 


CTD 


0.09 


184 


8394216-5181 


ONC 


Y 


Y 


N 


WTD 


0.36 


185 


8394454-5823 


ONC 


Y 


N 


Y 


WTD 


0.27 


188 


8396035-5861 


ONC 


Y 


N 


N 


WTD 


0.35 


198 


8399269-6612 


L1641 


N 


N 


Y 


WTD 


0.51 d 


203 


8406280-6293 


L1641 


N 


N 


Y 


PTD 


0.64 


204 


8406609-6247 


L1641 


Y 


N 


N 


CTD 


0.08 


211 


8411618-6426 


L1641 


Y 


N 


N 


CTD 


0.12 


218 


8416389-6503 


L1641 


Y 


N 


N 


CTD 


0.10 


219 


8416834-6225 


L1641 


Y 


Y 


N 


CTD 


0.00 


221 


8416978-6185 


L1641 


N 


N 


Y 


PTD 


1.18 


223 


8424570-6484 


L1641 


Y 


Y 


N 


CTD 


0.04 


227 


8440418-7404 


L1641 


N 


N 


Y 


PTD 


0.70 


229 


8444779-6608 


L1641 


Y 


Y 


N 


CTD 


0.04 
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Table 5 — Continued 



Trt 


name 


region 


n-K—6 


^13 — 31 


Hi VV yLufllll j 


subtype 










(A\ 

\V 


V>) 


w 


n\ 

V) 


\°> 


230 


8445558-6860 


L1641 


Y 


Y 


Y 


CTD 


-0.13 


237 


8457266-7161 


L1641 


Y 


Y 


Y 


CTD 


0.00 


239 


8464666-7838 


L1641 


Y 


N 


Y 


WTD 


0.36 


257 


8502614-7795 


L1641 


Y 


Y 


N 


CTD 


0.09 


260 


8508501-7431 


L1641 


Y 


Y 


N 


CTD 


0.02 


270 


8519428-8120 


L1641 


Y 


Y 


N 


CTD 


-0.25 


271 


8520557-7775 


L1641 


Y 


Y 


N 


CTD 


0.08 


275 


8535781-7830 


L1641 


Y 


Y 


N 


WTD 


0.37 


279 


8538911-7999 


L1641 


N 


Y 


N 


WTD 


0.31 


283 


8547523-7831 


L1641 


Y 


Y 


N 


CTD 


0.03 


290 


8562845-8151 


L1641 


N 


N 


Y 


PTD 


0.78 


291 


8564835-8250 


L1641 


Y 


Y 


N 


CTD 


0.04 


294 


8567670-8803 


L1641 


Y 


Y 


N 


CTD 


0.14 


297 


8570949-8577 


L1641 


N 


Y 


Y 


PTD 


1.14 


301 


8576836-8303 


L1641 


N 


N 


Y 


PTD 


0.63 f 


302 


8579480-8565 


L1641 


N 


Y 


N 


WTD 


0.39 


303 


8580634-8516 


L1641 


Y 


Y 


N 


CTD 


0.19 




04202606 


Tau 


N 


N 


Y 


PTD 


0.52* 




04125+2902 


Tau 


Y 


Y 


N 


CTD 


0.04 




CoKuTau 4 


Tau 


Y 


Y 


Y 


CTD 


0.06 




DMTau 


Tau 


Y 


Y 


Y 


CTD 


0.15 




GKTau 


Tau 


N 


N 


Y 


PTD 


0.99 




GMAur 


Tau 


Y 


Y 


Y 


WTD 


0.37 




HKTau 


Tau 


N 


Y 


N 


PTD 


0.77 




IPTau 


Tau 


N 


N 


Y 


PTD 


0.55 




LkCal5 


Tau 


N 


Y 


Y 


PTD 


0.56 




MH03 


Tau 


N 


N 


Y 


PTD 


1.10 




RYTau 


Tau 


N 


N 


Y 


PTD 


1.14 




UXTauA 


Tau 


N 


Y 


N 


PTD 


0.74 




V410 Xray6 


Tau 


Y 


Y 


Y 


CTD 


-0.04 




2MJ1124118 


Cha I 


Y 


Y 


Y 


CTD 


-0.50 




Baud43 


Cha I 


N 


N 


Y 


PTD 


0.45 c3 


T54 


CHX22 


Cha I 


Y 


Y 


N 


CTD 


-0.03 




CHXR22E 


Cha I 


Y 


Y 


N 


CTD 


0.18 




CRCha 


Cha I 


N 


N 


Y 


PTD 


0.71 


Til 


CSCha 


Cha I 


Y 


Y 


N 


CTD 


0.15 




IS091 


Cha I 


Y 


N 


Y 


CTD 


-0.22 


T25 


Szl8 


Cha I 


Y 


Y 


N 


CTD 


0.06 


T35 


Sz27 


Cha I 


Y 


Y 


N 


PTD 


0.44 c2 


T56 


Sz45 


Cha I 


Y 


Y 


Y 


WTD 


0.35 


T6 


SZCha 


Cha I 


Y 


Y 


N 


PTD 


0.61 




16126-2235 


Oph 


N 


Y 


Y 


PTD 


0.49 e3 




16201-2410 


Oph 


N 


Y 


Y 


PTD 


0.54 




16220-2452 


Oph 


N 


N 


Y 


PTD 


1.10 




16225-2607 


Oph 


N 


N 


Y 


PTD 


0.47 c2 



Table 5 — Continued 



ID 


name 


region 


nx-6 


"13-31 


EW(10(mi) 


subtype 


IDEF 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 






vjpn 


NT 
IN 




V" 
1 


r 1 U 


U.oo 




DOAR28 


Oph 


Y 


Y 


N 


CTD 


0.07 




GY195 


Oph 


N 


N 


Y 


PTD 


0.65 f 




IRS3 


Oph 


N 


N 


Y 


PTD 


1.38 




Rox44 


Oph 


N 


Y 


Y 


PTD 


0.90 




SR9 


Oph 


N 


N 


Y 


PTD 


0.53 


spID 50 


LAL 106 


N1333 


N 


N 


Y 


PTD 


0.80 f 


spID 52 


LAL 110 


N1333 


N 


N 


Y 


PTD 


0.55 


spID 57 


LAL 129 


N1333 


N 


N 


Y 


PTD 


0.57 


spID 73 


LAL 171 


N1333 


N 


N 


Y 


PTD 


0.72 


spID 101 


LAL 245 


N1333 


N 


Y 


N 


PTD 


0.75 f 


spID 110 


LAL 279 


N1333 


N 


Y 


N 


PTD 


0.81 


spID 116 


LAL 300 


N1333 


N 


N 


Y 


PTD 


1.71 


spID 136 


LAL 321 


N1333 


Y 


N 


N 


CTD 


0.03 


spID 137 


LAL 331 


N1333 


Y 


Y 


N 


CTD 


0.22 



Note. — Column (4) indicates whether an object passes the criteria of n^-_g; 
Column (5) indicates whether an object passes the criteria of ni3_3i; Column 

(6) indicates whether an object passes the criteria of EW (10/xm); Column 

(7) : subtype of TDs based on the values of IDEF in column (8). Sec Section 
4.1.1 for a detailed explanation of IDEF and the subtype separation points 
of IDEF.] Column (8): IDEF is the inner disk excess fraction. 

a OriA-44 and OriA-176 do not pass n#-_ 6 , 7113-31, and £W(10//m), but 
their SEDs resemble a TD's SED. 

b we did not calculate its spectral indices because its IRS spectrum is not 
complete, but its SED resemble a TD's SED. 

c OriA-8 and OriA-172 are selected as CTD, but IDEF > 0.25. The spectral 
type of the objects is M3 or later. When measuring ef using the median 
spectrum of objects with M3-M5 in Tau, the IDEF values are 0.45 and 0.65, 
respectively. 

d OriA-23 and OriA-198 are selected as WTD with IDEF = 0.51. The 
spectral type of objects is M0 and K6.5 respectively. 

cl OriA-149 and OriA-164 have IDEF<0.5, but they are selected as PTD 
because we do not know their spectral types, but the fluxes at wavelength 
ranges of 2-6 ^mare comparable to that of the median spectrum. 

c2 T35 (Sz 27) and 16126-2607 are PTDs, but their IDEF are slightly less 
than 0.5. Their spectral types are M0 and K7, respectively. 

c3 Baud 43 and 16126-2235 are PTDs , but their IDEF are less than 0.5. 
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Their spectral types are M3.25 and M3.0, respectively. When measuring IDEF 
using the median spectrum of objects with M3-M5 in Tau, the IDEF values 
are 0.72 and 0.71, respectively. 

f OriA-5, OriA-301, 04202606, GY 195, LAL 106, and LAL 245 have spectral 
types of M3 or later. Their IDEF estimated with K5-M2 median spectrum 
are higher than 0.5 and they are classified as PTDs. The IDEF estimated 
with M3-M5 median spectrum are much higher than the IDEF listed in this 
table: 0.94, 1.07, 0.87, 1.0, 1.5, and 1.14, respectively. 
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Table 6. number of transitional disks and age trend 



region 


Age (Myr) 


CTD 


WTD 


PTD 


total TD 


disks 


NGC 1333 


< 1 


2 





7 


9 


53 a 


ONC 


< 1 


15 


10 


7 


32 


127 b 


L1641 


1 


19 


5 


6 


30 


114 b 


Tau 


1.5 


1 


1 


8 


13 


154 c 


Chal 


2 


6 


1 


1 


11 


68 d 


Oph 


0.3-2.1 


1 





9 


10 


90 e 



Note. — Number of transitional disks and the median age of each 
star-forming region. The median age of each star-forming region is taken 
from the literature listed below and references therein. The last column: 
n(disks) is the number of objects in our sample, which is complete for Class 
II objects and transitional disks with host stars of type M4 or earlier (see 
§2). 

a We keep 53 objects after taking 4 Cla ss I objects out from 
the 57 disk dominant objects reported in lArnold et al. ( 2012 ). 



b from this work. 

c This number c omes from IRS survey of Taurus by 
Furlan et all |201ll) . 



d This number co mes from IRS survey of Chamaeleon I by 
Manoi et all (|201lh . 



°This number is the number of disk dominant objec t s iden ti- 
fied from IRS survey of Ophiuchus bv lMcClure et al.l (|2010l ). 
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Table 7. Correlation summary I 



Correlation 


IN 


OL 


o 
P 


COTT° 


f (/o) 


p d 

Grange \y0) 


M vs. M,/M 


69 


-8.5±0.1 


1.6±0.3 


0.5±0.1 


0.001 


<0.001-0.1 


L x vs. M*/M 


56 


30.1±0.1 


1.1±0.3 


0.5±0.1 


0.01 


<0.001-0.2 




55 


30.1±0.1 


1.1±0.3 


0.4±0.1 


0.3 


0.01-2.6 


M vs. L x 


42 


-14.6±7.7 


0.2±0.3 


0.2±0.2 


20.4 


0.9-100 


M*/M Q VS. Rwall 


98 a 


-l.OiO.l 


0.7±0.1 


0.8±0.1 


<0.001 


<0.001 




97 


-l.OiO.l 


0.7±0.1 


0.7±0.1 


<0.001 


<0.001 


M VS. 


69 


-10.1±0.4 


1.1±0.4 


0.4±0.1 


0.07 


0.001-1.2 


L X VS. R^aii 


56 a 


28.5±0.3 


1.2±0.3 


0.6±0.1 


<0.001 


<0.001-0.01 




55 


28.5±0.4 


1.2±0.3 


0.5±0.1 


0.01 


<0.001-0.3 


M/M[M±,R W ] vs. R wa H 


69 


0.05±0.3 


-0.7±0.3 


-0.3±0.1 


1.2 


0.1-10 


L X /L X [M*,R W ] vs. i? waii 


56 a 


0.05±0.3 


-0.2±0.3 


-0.1±0.2 


46 


2.5-100 




55 


0.1±0.4 


-0.2±0.3 


-0.1±0.2 


47 


2.6-100 


M/M[L X ] vs. L x 


42 


28.6±5.9 


-1.0± 0.2 


-0.7±0.1 


<0.001 


<0.01-<0.001 


L* vs. M-k 


98 a 


0.3±0.1 


1.6±0.1 


0.8±0.04 


<0.001 


<0.001 




97 


0.3±0.1 


1.5±0.1 


0.8±0.04 


<0.001 


<0.001 


L*/L Q VS. 


98 a 


-1.9±0.1 


1.6±0.1 


0.9±0.04 


<0.001 


-C0.001 




97 


-1.8±0.2 


1.5±0.1 


0.9±0.04 


<0.001 


<0.001 


Lx vs. L*/L 


56 a 


29.8±0.1 


0.8±0.2 


0.6±0.1 


<0.001 


<0.001-<0.01 




55 


29.9 ± 0.1 


0.8±0.2 


0.6±0.1 


<0.001 


<0.001-0.01 


Lx/L* vs. M*/Af© 


56 a 


-3.8±0.1 


-0.4±0.4 


-0.2±0.1 


14 


2.5-46 




55 


-3.8±0.1 


-0.3±0.3 


-0.2±0.1 


14 


2.6-47 


M vs. L*/L Q 


69 


-8.8±0.1 


0.8±0.2 


0.5±0.1 


<0.001 


<0.001-0.1 


M vs. Lx /L* 


42 


-10.2±1.1 


-0.4±0.3 


-0.3±0.2 


5.4 


0.1-53 



Note. — Summary of the results from the Bayesian linear regression between X and Y: Y vs. X; 
logY = a + filogX. Objects has a known spectral type are only included for the correlations. 

a The results when we include OriA-88, which has a large value of R wa ii, M±, and L*, 
and therefore is an outlier. 

b corr is the linear correlation coefficient between logY and logX. 

C P is the probability of getting corr in random distribution. 

d Prange is the range of P to get corr in la uncertainty. 



Table 8. Correlation summary II: details by subtypes of TDs 



Correlation 


TD type 


N 


a 


/3 


corr h 


P c (%) 


P d /0/\ 


MRP 


XPE f 


GPFs 


M vs. Nh/M Q 


CTD 


34 


-8.9+0.6 


0.7+1.8 


0.1+0.3 


57 


1.9-26 






Y 




WTD+PTD 


35 


-8.2 ± 0.1 


1.9+0.5 


0.7+0.1 


<0.001 


<0.001-0.01 






Y 


L x vs. M*/M Q 


CTD a 


30 


30.3+0.2 


1.8+0.6 


0.6+0.2 


0.05 


<0.001-2.9 






Y 




CTD 


29 


30.5+0.3 


2.2+0.9 


0.5+0.2 


0.6 


0.002-11 






Y 




WTD+PTD 


26 


30.0+0.1 


1.0+0.4 


0.5+0.2 


0.9 


0.007-14 






Y 


M vs. R waU 


CTD 


34 


-9.1+0.6 


-0.1+0.6 


-0.0+0.2 


100 


26-100 


N 


N? 


Y 




WTD+PTD 


35 


-10.8+0.6 


1.9+0.5 


0.7+0.1 


<0.001 


<0.001-0.01 


Y 


Y? 


Y 


M/M(AL) vs. R wall 


CTD 


34 


0.8+0.7 


-1.7+0.6 


-0.6+0.2 


0.02 


<0.001-1.9 


N 


Y? 


Y 




WTD+PTD 


35 


-0.9+0.6 


0.2+0.5 


0.1+0.2 


57 


8-100 


N 


N 


Y 


L x vs. RwaU 


CTD a 


30 


28.3+0.4 


1.4+0.3 


0.7+0.1 


0.002 


<0.001-0.05 




Y? 


Y 




CTD 


29 


28.2+0.4 


1.5+0.4 


0.7+0.2 


0.002 


<0.001-0.6 




Y? 


Y 




WTD+PTD 


26 


28.9+0.9 


0.7+0.7 


0.3+0.3 


14 


0.1-100 




Y? 


Y 


L X /L X {M+) vs. R wall 


CTD a 


30 


-0.2+0.3 


0.1+0.3 


0.1+0.2 


60 


11-100 


N 


N 


Y 




CTD 


29 


-0.2+0.4 


0.1+0.4 


0.1+0.2 


61 


11-100 


N 


N 


Y 




WTD+PTD 


26 


0.7+0.9 


-0.8+0.7 


-0.3+0.3 


14 


0.1-100 


N 


N 


Y 


M vs. L x 


CTD 


24 


-17.1+9.9 


0.3+0.3 


0.3+0.3 


15 


0.2-100 


N 


N? 


Y 




WTD+PTD 


18 


-8.3+10.2 


0.0+0.3 


0.0+0.3 


100 


23-100 


N 




Y 


M/M{L X ) vs. L x 


CTD 


21 


28.1+10.0 


-1.0+0.3 


-0.6+0.2 


0.2 


<0.001-5.3 


N 


Y? 


Y 




WTD+PTD 


18 


35.7+9.7 


-1.2+0.3 


-0.7+0.2 


0.1 


<0.001-3.5 


N 


Y? 


Y 


M*/M vs. Rwi 


CTD 


46 


-1.0+0.1 


0.6+0.1 


0.8+0.1 


<0.001 


<0.001 




Y 






WTD+PTD 


52 


-1.1+0.2 


0.8+0.1 


0.7+0.1 


<0.001 


<0.001 




Y 




A/ vs. L*/Lq 


CTD 


34 


-9.3+0.2 


-0.3+0.4 


-0.2+0.2 


26 


1.9-100 






Y 




WTD+PTD 


35 


-8.6 ± 0.1 


1.2+0.2 


0.8+0.1 


<0.001 


<0.001 






Y 


vs. L*/L e 


CTD a 


30 


29.9+0.1 


0.8+0.2 


0.6+0.2 


0.05 


<0. 001-2. 8 


Y 




Y 




CTD 


29 


30.0+0.1 


1.0+0.3 


0.6+0.2 


0.06 


<0. 001-3. 2 


Y 




Y 




WTD+PTD 


26 


29.7+0.1 


0.8+0.3 


0.6+0.2 


0.2 


<0.001-4.3 


Y 




Y 



Note. — Summary of the results from the Bayesian linear regression between X and Y: Y vs. X; logY = a + filogX. 
The index a > b > c > and d have the same meanings as in Table [7] 

C MRI: The inside-out disk evolution driven by Magnetic Rotational Instability. 
l XPE: The disk dispersion by photoevaporation, specially X-ray photoevaporation 



e GPF: Dynamical effects by infant giant planet formation. 

Note. — N or Y checks indicate if a mechanism in the last three columns can explain the corresponding trends in a row. 



I 
I 
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Table 9. Correlation summary III: TDs in Orion A 



Correlation 


N 


a 




corr b 


P c (%) 


Prange (%) 


M vs. Mi,/M @ 


50 


-8.6±0.2 


1.3±0.5 


0.4±0.2 


0.4 


<0. 001-16 


L x vs. M*/M 


35 a 


30.2±0.1 


1.0±0.3 


0.5±0.2 


0.2 


<0. 001-8 




34 


30.2±0.2 


0.9±0.4 


0.4±0.2 


1.9 


0.02-26 


M vs. L x 


33 


-20.3±8.7 


0.4±0.3 


0.3±0.2 


9 


0.3-58 


M*/M© vs. R waU 


55 a 


-l.liO.l 


0.7±0.1 


0.8±0.1 


<0.001 


<0.001 




54 


-1.0±0.1 


0.6±0.1 


0.8±0.1 


<0.001 


<0.001 


M vs. R wa u 


50 


-9.8±0.5 


0.7±0.4 


0.3±0.2 


3.4 


0.2-49 


L x vs. R wa u 


35 a 


28.8±0.3 


1.0±0.3 


0.6±0.1 


0.01 


<0.001-0.2 




34 


28.8±0.4 


1.0±0.4 


0.5±0.2 


0.3 


<0.001-8.5 


M/M\M*\ vs. R wall 


50 


0.3±0.3 


-1.0±0.3 


-0.5±0.1 


0.02 


<0.001-0.4 


Lx/LxlM^Rw] vs. R wall 


35 a 


0.3±0.3 


-0.3±0.3 


-0.2±0.2 


25 


1.7-100 




34 


0.4±0.4 


-0.3±0.4 


-0.2±0.2 


26 


1.9-100 


M/M[L X ] vs. L x 


33 


24.0±6.1 


-0.8±0.2 


-0.7±0.1 


<0.001 


<0.001-0.02 


L x vs. 


35 a 


30.0±0.1 


0.6±0.2 


0.6±0.1 


0.01 


<0. 001-0. 2 




34 


30.0 ± 0.1 


0.7±0.2 


0.6±0.2 


1.0 


<0.001-1.9 


Lx/L* vs. M*/M© 


35 a 


-3.8±0.1 


-0.7±0.3 


-0.4±0.2 


1.7 


0.01-25 




34 


-3.8±0.1 


-0.6±0.4 


-0.3±0.2 


8.5 


0.3-57 


M vs. L,/L Q 


50 


-8.9±0.1 


0.6±0.3 


0.4±0.2 


0.4 


<0. 001-16 


M vs. L x /L* 


33 


-8.3±1.3 


0.1±0.4 


0.1±0.2 


58 


9-100 



Note. — Summary of the results from the Bayesian linear regression between X and Y: Y vs. X; 
logY = a + /3logX. Objects has a known spectral type and in Orion A are only included for the 
correlations. 

The index a > b > c i and d have the same meanings in Table]?] 



Table 10. Correlation summary IV: details by subtypes of TDs in Orion A 



Correlation 


TD type 


N 


a 


P 


corr h 


P c (%) 


Pran e d (%) 

range 


M vs. M*/M Q 


CTD 


30 


-9.2+0.5 


-0.2+1.3 


0.0+0.3 


100 


11-100 




WTD+PTD 


20 


-8.3 ± 0.3 


2.1+0.7 


0.6+0.2 


0.5 


0.002-8.1 


L x vs. M+/M Q 


CTD a 


24 


30.3+0.2 


1.4+0.5 


0.5+0.2 


1.3 


0.01-15 




CTD 


23 


30.5+0.4 


1.9+1.2 


0.4+0.2 


5.7 


0.2-36 




WTD+PTD 


11 


30.2+0.1 


0.5+0.4 


0.5+0.3 


12 


0.3-56 


M vs. R waU 


CTD 


30 


-8.8+0.6 


-0.3+0.6 


-0.1+0.2 


60 


11-100 




WTD+PTD 


20 


-11.1+0.9 


2.0+0.8 


0.7+0.2 


0.06 


<0.001-2.5 


M/M(M+) vs. R wall 


CTD 


.30 


1.0+0.6 


-2.0+0.6 


-0.7+0.2 


0.002 


<0.001-0.5 




WTD+PTD 


20 


-1.1+0.8 


0.2+0.7 


0.1+0.3 


67 


8.1-100 


L x vs. R wa n 


CTD a 


24 


28.6+0.4 


1.1+0.4 


0.6+0.2 


0.2 


•CO.001-5.3 




CTD 


23 


28.6+0.6 


1.2+0.5 


0.5+0.2 


1.5 


0.02-16 




WTD+PTD 


11 


29.5+0.6 


0.5+0.5 


0.4+0.4 


22 


0.3-100 


L x /L x (hU) vs. R wall 


CTD a 


24 


0.2+0.4 


-0.2+0.30 


-0.1+0.2 


64 


15-100 




CTD 


23 


0.2+0.6 


-0.2+0.6 


-0.1+0.3 


65 


5.9-100 




WTD+PTD 


11 


1.1+0.6 


-0.8+0.5 


-0.6+0.3 


5.1 


0.02-37 


M vs. L x 


CTD 


23 


-14.3+10.5 


0.2+0.4 


0.1+0.3 


65 


5.9-100 




WTD+PTD 


10 


-37.3+29.7 


1.0+1.0 


0.5+0.4 


14 


0.04-78 


M/M(L X ) vs. L x 


CTD 


23 


31.0+9.0 


-1.1+0.3 


-0.7+0.2 


0.02 


<0.001-1.5 




WTD+PTD 


10 


5.8+22.9 


-0.2+0.8 


-0.2+0.4 


58 


6.7-100 


M,/M vs. R wall 


CTD a 


33 


-1.0+0.1 


0.6+0.1 


0.8+0.1 


<0.001 


<0.001 




CTD 


.32 


-0.9+0.1 


0.5+0.1 


0.7+0.1 


<0.001 


<0.001-0.03 




WTD+PTD 


22 


-1.1+0.2 


0.8+0.2 


0.8+0.1 


<0.001 


<0.001-0.03 


M vs. L*/Lq 


CTD 


.30 


-9.3+0.2 


-0.4+0.4 


-0.3+0.2 


11 


0.5-60 




WTD+PTD 


20 


-8.8 ± 0.2 


1.4+0.4 


0.8+0.1 


0.002 


<0.001-0.06 


L x vs. L*/Lq 


CTD a 


21 


30.1+0.1 


0.7+0.2 


0.6+0.2 


0.2 


<0.001-5.3 




CTD 


23 


30.0+0.1 


0.7+0.2 


0.5+0.2 


1.5 


0.02-16 




WTD+PTD 


11 


30.0+0.1 


0.5+0.2 


0.6+0.3 


5.1 


0.02-37 



Note. — Summary of the results from the Bayesian linear regression between X and Y: Y vs. X; logY = 
a + filogX. Objects has a known spectral type and in Orion A are only included for the correlations. 
The index a ' f, > c > and d have the same meanings in Table [7] 



